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A compact and light gear box is one of the utmost 
consideration in the design of machinery space on board a 
Naval Ship. Compactness is achieved by a rational and systematic 
design of the gearing so that it occupies the minimum possible 
floor area. A mathematical model for the compactness of double 
reduction gearing set is formulated. Interior Penalty Function 
approach is used to arrive at the optimum size of the gears. 

The effectiveness of the proposed algorithm Is proved by 
comparing the computer-aided design with the existing one. 
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CHAPTER 1 
INTRODUCTION 

Not long ago, weight was a critical factor mainly 
in Air Craft design. Today, weight reduction is the top 
design goal for a wide range of products starting from the 
24' Marine gears to small wrist watches. Weight reduc- 
tion usually means volume reduction, which in turn lowers 
cost of materials, handling and labour charges. 

The designer of a Marine reduction gear box is 
mainly concerned with the above problem. The lighest gear 
box occupying the minimum deck space is considered the best 
design* The marine gears are required to take heavy and 
sudden loads, as in time of emergency, the ship is required 
to shoot to the maximum speed (say 25 knots) from stand still 
in an interval of 20-30 seconds thereby causing excessive 
impact on the gear teeth. The second requirement for 
naval gears is the quiet operation of the gear box to 
avoid detection of the ship’s presence, due to the water 
borne noises, by the underwater operating submarines and 
other sophisticated gadgets employed by the enemy. And at 
the same time it is desirable to operate one’s own ship 
quietly enough, to be able to hear water borne noises of 
the enemy's ships. 
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The bottom of a ship is generally of girder and 
steel plate construction and is therefore relatively • 
flexible - a condition which may affect the design of the 
gear box. The position of the propeller shaft may influence 
the design by limiting the size of the bull wheel that can 
be accommodated between the shaft centre and the ship l s 
bottom. 

The first limitation in the design is met by 
providing flexible shafts and second by having the speed 
reduction in 2 stages, as the bull wheel in the 2 stage 
reduction is much smaller than, the one in single reduction. 
Quietness in operation is achieved by making the gears with 
large number of teeth and cutting the teeth with extreme 
accuracy. The problem of marine gear design is complete 
if the above requirements are met and the components are 
designed for adequate power transmission* 

1.1 Critical Review of Past Investigations 

Ray C. Johnson^ , Staff Engineer, IBM, New York, 
has put forward a procedure for the optimum design of heli- 
cal gears. But the author himself is not very sure of the 
reliability and accuracy of his method. His design procedure 
has certain inherent drawbacks. In his procedure of 
optimum design, he has incorporated all significant 
information in the equation form, and much ox this 
information was originally available only in graphical or 
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tabular form, suitable for analysis or a cut -and -try type 

of design procedure. Secondly the author had put certain 

range restrictions on the parameters, constituting his 

design equations, which restrict the scope of the proposed 

procedure to be used for optimal design, 

( 20 ) 

VJ.ik. Tuplin had suggested some curves for 
determining helix angles for minimum centre distance 
between the gear sets for various gear ratios and shaft 
angles. The method suffers from the defect that one has 
to decide the module and the gear ratios in the inter- 
mediate stages arbitrarily. The design so obtained can 
be far away from the optimal design as the designer 
would not know what specific combination of module and 
gear ratio will yield the optimal gear set. 

E.J. Weliauer^ had developed an analytical 
procedure using a number of strength rating formulae 
based on the wear rating formula developed by AGMA 
standards. The author claims that the design procedure 


obtained from his equations and charts result in an 
optimum design from the consideration of space occupied 
by the gear set. Though his method is workable for- 
ge ars transmitting moderate horse powers but suffers 
from almost the same drawbacks as Johnson’s optimal 
design procedure. 
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( i9) 

R.J . Wills has suggested some new equations 
and charts which pick up the lightest weight spur gears. 

The method presented indicates gear sizes and ratios that 
will permit lightest possible design while still meeting 
the horse power and surface durability requirements* 

Though the equations and charts can be extended for 
designing helical gears but the method cannot incorporate 
the special requirements of marine gearing design and is 
suitable only for conventional gear design* 

C-.J. Huebne r v has developed a computer 
oriented procedure for the design of gears. Balanced 
life concept of gear design is outlined in which gear 
and pinion are designed to fail simultaneously. The 
author has solved an example shox^ing how the balanced 
life concept allows combination of minimum size and 
maximum capacity. But again the method suffers from the 
fact that special design requirements cannot be incorporated 
in this method. Moreover it cannot be used for multistage 
reduction gearing as the author does not suggest any 
criteria to choose the gear ratios in the intermediate 
stages of reduction. 

1*2 Present Investigation - Its Scope 

In the present dissertation a mathematical 
model of two stage marine reduction gear box has been 
developed. The algorithm takes into consideration the 
various features of marine gearing. All the component 
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gears have been designed for the maximum Horse Power, 
For this, a computer programme has been developed. 
Important correction data for the gear teeth has also 
been generated in the same programme , 

In view of the fact that the compactness of the 
gear box loads to economy of space, machining time and 
weight, the minimum area occupied by the gear box has 
been chosen as the optimality criteria. 

The optimal design problem is formulated as a 
mathematical programming problem. The optimisation 
algorithm chosen is the Fletcher and Powell's Variable 
Metric method. The method is iterative in nature and 
seeks the optimum from within the constraint set. To 
prove the effectiveness of the proposed algorithm, 
three example problems were taken up. First one refers 
to the reduction gearing for a Naval destroyer snip 
(40,000 H.P.), the second one is for a cargo ship 
(60,000 H.P.). In the third problem, tic gears have 
been designed to transmit 75000 H.P., over a speed 
reduction ratio of 20.0. 
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CHAPTER II 

MATHEMATICAL FORMULATIOH 

In this Chapter, the design variables used in the 
mathematical modelling of the gear set, are introduced 
and the objective function. is formulated. The design 
constraints are expressed in the form of inequalities and 
a penalty 'function to account for the constraints is 
generated. Also the solution scheme to solve the problem 
is discussed in brief. 

2.1 Design Variables 

The numerical quantities for which values are to 
be chosen In producing a design will be called 'Design 
Variables'. Other names have also been used, and, of 
course, the terminology used here is merely the author's 
preference. Perhaps the most descriptive of the alterna- 
tive names are 'Construction parameters' or ^Construction 
Variables' . 

The problem uses 11 variables. The design 
variables are the number of teeth on pinion in both the 
stages, the modules for both the stages, the transmission 
ratio in the first stage (the transmission ratio in the 
second stage being fixed from the consideration of overall 
transmission ratio ) 5 the helix angles for both the stages, 
the face width coefficients, the angle ©( (used for positi- 
oning the gear box centrally), and the error' term which 
is the allowable manufacturing error for the gears. 
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It has been tactly decided to reduce the turbine 
speed to the propeller speed in two stages only. Experience 
■in the field of marine gears confirms that it is advisable 
to reduce the speed in two stages only for a speed reduc- 
tion ratio of the order of 15-30. 

The pressure angle 0 for the gears have been 
arbitrarily taken as 20 degrees. Here again it has been 
found that most of the heavy duty marine gears use 20.0° 
pressure angle because of the practical considerations. 


2.2 Nomenclature of Important Components 

The gears are named in single subscript. Zp-j_ 


and Zpg stand for the number of teeth on pinion in the 
first and the second stage respectively. Similarly Zg^, 
Zg 2 > SIA1, SIA2, C rfl , stand for the number of teeth 

on gears, the helix angles, the face width coefficients 
for the first and the second stages respectively. FM0B1 
and PM0D2 stand for the modulo s in the same way. 


2.3 Design Space 

If the number of independent variables in a 
particular configuration of a gear box is n, then the 
n-dimensional Eucledian space is known as the design space. 
Every point in the space represents a design even if it 
is absurd (as negative number of teeth or a gear digging 
into the other shaft). The movement from point to point 
in this space is represented by a vector 
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2.4 Constraint Surface 

The set or locus of values of X which satisfy 
tho equation g^(X) = 0 form a surface in t he design space. 
This is not a 2-dimensional sub space, but an n-1 dimensional 
on’o where n is the number of design variables. It is a 
surface in the sense that it outs the space into two 
regions, one where g 0 (the acceptable region) and the 
other where g 4 0 (the unacceptable region). Hence the 
surface that separates these two regions is called a 
constraint surface. 

The portions of the respective constraint surfaces 
which bound the feasible region form a patch work of the 
constraints called the composite constraint surface* 

Points within this region (i.e. where g^(X) j> 0), are 
called free points or unconstrained designs and points 
on the surface (i.e* designs which are such that at least 
one g^(X)=0) are called bound points or constrained designs. 
It is possible for the acceptable regions or tho feasible 
region to be composed of two or more disjoint subregions. 

Tho sub-space where two or more gj_(X)=0 will be 
referred to as an intersection. die dimension of an 
intersection is n-r where r distinct constraints intersect. 

2.5 Objective Function 

Of all designs which arc feasible some will be 
more desirable or better than others. If this is true then 
there must be some quality which the better designs have 
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more of, than do the less desirables ones have. If this 
quality can be expressed as a computable function of the 
design variables 5 we can consider optimizing to obtain a 
"Best Design” ♦ The function, with respect to which the 
design is optimized, is called the objective function# We 
will designate the objective function as F(X) . The best 
feasible design is one for which 

g±(X) ^ 0 , i = 1,2,3, -.N c 

and 

F(X) is optimum. 

The selection of an objective function for an 
engineering design problem can be one of the most important 
decisions in the whole optimum process. In some situations, 
an obvious objective function exists because the need which 
the design is to fulfil is better served by some designs 
than others. In some design situations, there may appear 
to be two or more quantities which should be objective 
functions. Such situations are usually handled in one 
of the three ways ; 

(1) Composite objective function is formulated 

(2) A limit is set for one of the functions and it is 
used as a constraint 

(3) A trade off study is done* 

In the present work the area occupied by the 
triangle formed by joining the centres of the three shafts 
is chosen as the optimality criterion (see fig. 2 . 1 ). 
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The following are the steps used to calculate 
the objective function ; 


TRAN2 

ss 

20.0/TRAN 1 

D P1 

3 

V /(I W Cos l f ; i ) 

D p2 

zz 

Z p2 /(P nd2* Cos , Y 2 ) 

b 

1 

rz 

C nfl* 7X * Cos ■■■ y p ndl 

b 2 

3 

C nf2* 7V * Cos X^ 2 /P nd2 

V 

= 

D pl . THAN 1 

V 


D p2 . TRAN 2 

CD 1 

= 

D pl . (TRAN 1 + 1.0)/2 .0 

CD 2 

r 

D p 2 . (TRAN 2 + 1.0) /2.0 

CD33 

= 

GDI 2 + CD2 2 - 2CD1 . CD2 

CD3 

Cos S 

W 3 

\/CD3S 

CDl 2 + CD3 2 - CD2 2 


2CD1 . CD 3 

(3 


AH Cos (Cos g ) 

FO 


CDS . CDl . Sin (5 


2.0 


2,6 Normalised Constraints 

The general form of expressions of design constraints 
is 

g x (X) = 0.01 - Fl(X) ^ 0 

g 2 (X) = F2(X)- 10,000 > 0 

where it is required that 
FI (X) ^0.01 
F2 (X) ^10,000 
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The above is an example where the magnitude and 
sensitivity of the two constraints with respect to change 
in design variables are quite different* This results in 
difficulty to arrive at the optimum point efficiently. The 
disparity in magnitudes and sensitivities of the design 
constraints can be avoided by normalising them, i.e., by 
forcing them to take values between 0 and 1. If the 
inequality constraint expresses the difference between 
two variables, then it should be normalised with respect 
to the variable of greater value. Sven when the magnitudes 
of the constraints are controlled, their sensitivities 
may vary widely. Nevertheless, normalising always improves 
the handling of constraints. 

2.7 Constraints 

The design restrictions, the satisfaction of 
which distinguishes acceptable designs from the unaccep- 
table designs, will be called collectively constraints. 
There are two useful categories of constraints in engineer- 
ing problems. This grouping is not necessarily definitive, 
as it may not always bo easy to classify constraints along 
these lines. 

A constraint which restricts the range of design 
variables for reasons other than the direct consideration 
of performance of the design will be called a side cons- 
traint. Constraints which derive from those performance 
or behaviour requirements which are explicitly considered 
are called behaviour constraints. 
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Another type of constraint which arises in some 
engineering problems is that of the discrete valued design 
variables. In such cases the design variable is not to 
be selected from a continuous range of values but is 
permitted to take on only one of a discrete set of values* 
For example in the gear design problem one of such cons- 
traints is that the nodules can only take some predetermined 
preferred values which do not form a continuous set. 

The following is the list of side and behaviour 
constraints for the problem. 


SIDE CONSTRAINTS : The face width coefficients for the 
marine gears vary from 4.0 to 12.0, the following constraints 


take care of this : 

GX (5) « 


1.00 


C nfl 

’ TSTO 


GX (6) 
GX (7) 


1.00 

1.00 


4.00 

C nfl 

g nf2 

12.0 


GX (8) 


1.00 


4.00 

C nf2 


The minimum, and the maximum values recommended 
for the helix angles are 15.0° and 35.0° respectively. 
Hence the constraints : 


GX (9) 


1.00 



35.0 


GX (10) = 1.00 - 3*4^ 

\f>i 
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GX (11) = 1,00 - 
GX (12) = 1.00 - 

r 2 

To achieve quietness in operation, a very high 
degree of accuracy in the manufacture of gears is desirable# 
The following constraints impose restrictions on the 
allowable error in manufacturing. 


±2 

35,0 


GX (2l) = 1.0 - 


GX (22) 


1.0 - 


0,0005 
ERR . 

ERR 


0.0012 

The most commonly used values for the length of 
the face varies from 1,5 to 6.00 P^g. Therefore the 
constraints : 


GX (25) = 1.00 - 


GX (26) = 1.00 - 


GX (27) = 1.00 - 


GX (28) = 1.00 - 


6 *°* P ndl 
1.5 -P 


Ml 


"1 

^2 


6.0 *P. 


nd2 


1.5 *P. 


MS 


To avoid too high a reduction taking place in the 
first stage alone, the following constraint has been 


imposed : 


GX (35) = 1.0 - 


TRil'Tl 

7.00 
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BEHAVIOUR CONSTRAINTS : lo avoid interference of the gear 
teeth, the following constraints on the pinion teeth have 
been imposed t 

GX (1) =1.0 - 

Zpl 

GX (2) = i.O - ^ L 

Z p2 

The following constraints on the face contact 
ratio and the profile contact ratio have been imposed to 
assure positive control of the gear sets for all possible 
angular positions t 


GX (3) 


1 . 0 ~ 


FOR 1 


0.555 PCR1 
FOR 1 


GX (4) 


, 1 , 0.555 PCR 2 

1 '°lcR 2 + FOR 2 


The following constraints have been incorporated 


to make sure that the individual gears can transmit the 
designed H.P. in wear and bending. 

GX (13) = 1.0 - 


GX (14) = 1.0 - 


WD2 

WW2 


02 (15) = i* 0 - frn 


GX (16) = 1.0 - 


SE. 


HFtf2 


GX (17) = 1.0 - 


HP 

HPB1 
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The following constraints check the designed 
gears for adequate power transmission with the k.G.M.h, 
standards. ; 


GX (IS) =1.0 


licit £±UlCjc ll 

Sx J'B 1 « ^13_*^hl 


GX (20) = 1.0 


i?c2* G t 2* C r2^ 

s;iB2 c ? 2 .c h2 


GEOMETRY CONSTRAINTS s The following constraints prevent 
the gears from intersecting with the shafts of the 
adjoining gears s 


GX (23) = 1.0 - 


GX (29) = 1.0 - 


GX (30) = 1.0 - 


GX (31) = 1.0 - 

GX (32) = 1.0 - 


0.6 Dg2 

CD 3 

2.0 » DSH1 

V 

2.0 * DSH2 

D p2 

2.0 ♦ DSH3 


^p2 

2.0(CD1 -15.0) 


The following three constraints have been incorpo 
rated so that the gearing set is laid out as centrally 
as possible ! 

IBS (0.5 CD3 - CD1 • Cos$ ) 

GX (24) = 1.0 - 


0.25 CD3 



16 


GX (33) 


1.0 - 


90.0 


GX (34) = i.o *. 

2.7-1 Positivity of Design Variables 

Negative values for the number of teeth, trans- 
ission ratio ana the moaules of the gears are meaningless 
To restrain the minimisation algorithm from going into 

this region oi design space, positivity constraints are 
added. 


G,(X) = 


where 


h ~ 2Wl 

^max^i" ^rain^i 




r lj2 ...... ? I'T 


(2.3) 


N - number of independent variables 
X ± = i th design variable 
^min^i = minimum value of ith design, variable 
(x maz^i = maximum value of i th design variable. 


2.7-2 Gear Teeth 

The gear has to have a minimum number of teeth 
to avoid involute interference. The following inequality 
if satisfied assures that there will be no undercutting 
of gear teeth. 

4 m /„ 

— ■§— tN + m) ^ 2 W N + N 

Sm 0n & V ng np np 


2 
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where 


“ " 1S a constant (Addendum Constant) of propor- 
tionality depending upon the standard gear 
tooth system. It equals 0.8 for 20° stab tooth. 

0 n = pressure angle referred to the normal plane 
N np = number of teeth in the formative gear set 
for the pinion 

N ng = number of teeth in the formative gear set 
for the gear 

On simplification we get 

w 2 -*/ 4ra(aK np + m) 

np r — s- 

(2a+l) Sir. 0 n 

where 


we get 


a - Gear Ratio 

nssummg a to be any high value say infinity 


N np^' 


2 ( 0 , 8 ) 

(Sin 20°) 2 


13.7 


A-/ 


14.0 


This constraint is applied to all the gears in 
the system and is expressed as 


GX (1) 1.0 - 

Z P1 

(2.4) 

GX (2) * 1,0 - ^2 

2p2 
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2 .7-3 The Contact Ratio 

Contact ratio is an important consideration in 
kinematic studies of helical gears. There are two signi- 
ficant types of contact ratios used in the terminology of 
helical gearing. 


PROFILE CONTACT RATIO : The profile contact ratio - PCR, 
is defined as the arc of action divided by the circular 
pitch, both being measured in diametral plane. The profile 
contact ratio Is an index for the average number of teeth 
in contact in any fixed diametral plane. The PCR is a 
function of the pitch helix angle , the number of 
pinion teeth N p ; the gear ratio a, the normal pressure 

angle 0 n and the addendum. Expressed mathematically PCR 

C 
P 


, . ___ 0.032 „ 0.104 „ 1.504 

is given by FOR = 0.555 • a . N n . Cos ^ 


FACE CONTACT RATIO : It is defined as the ratio of helical 
advance of a tooth element across the active face width to 
the circular pitch, both being measured in diametral plane. 
Expressed mathematically, face contact ratio FCR is 
defined by 


b tan ‘‘I s 

FCR = ~ 

Pc 

b = Active face width 
P c = Circular pitch 

In order to assure positive control of the gear 
set for all po-ssible angular positions, a definite 
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relationship between the contact ratios must be satisfied. 
From a kinematic consideration of helical gear sets, the 
stipulation to assure positive control of the gear set 
for all possible angular positions is 

b tan ^ P c - arc of action 
b tan ^ arc of action 

Pc ^ i - P ~ 

°r FCR > i - for 

Preferably for good control and smooth action, 

FCR ^ i - for 

It is often entirely satisfactory to have, a 
profile contact ratio PCR less than unity provided FCR is 
large enough to satisfy the above equation. 

2.7-4 Wear Strength 

The load carrying capacity of helical gears is 
normally limited by pitting resistance of the tooth surface 

(wear). This resistance is commonly called the ’durability 
capacity. 1 

P itting - Pitting is defined as surface fatigue failure 
of the material caused by repeated surface or sub-surface 
stresses that exceed the endurance limit of the material. 

It is characterised by the removal of metal and the forma- 
tion of cavities. Pitting can be of two types : 

Sl4iial_pJiting - This may occur at the beginning 
ope ration and continue only until the over stressed 
localized high areas of the surface have been reduced, 


thus 
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providing sufficient contact area to carry the load 
Without further deterioration* Such pitting is not 
Serious, because it is corrective and nonprogressive » 

Ssi. truc t l.ve Pitting - This usually starts below 

the pitch line. The size and number of pits increase 
until Smoothness of operation is impaired* The remaining 
surface fails in a similar manner, and finally the tooth 
shape is destroyed. The pits constitute stress raisers 
which can lead to failure by fatigue breakage* 

Corrective and nonprogressive initial pitting 
is not considered serious. The AGMA formulae used here 
in the design are based on destructive pitting failures-. 

FAILURE - Pitting failure can be caused by any one or 
a combination of surface and subsurface shear, tensile 
and compressive stresses, depending on the specific 
distribution of inclusions, microstresses or other 
localized conditions. 

Since the calculated and allowable stresses have 
a linear relation^ the exact stress and location need not 
be determined. Instead the basic equation for contact 
stress is expressed as a 'contact stress number'," S c 
Its numerical value is equal to the surface compressive 
stress and. is 3-§- times the maximum subsurface shear 
stress. The contact stress number is 
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c m .c. 


) (“1 


) 


The relation "between calculated and allowable contact 
stress number is 


S c <^SALB 



) 


The .power capacity of gear set on wear basis is 


HPW = 


V i 

126,000 


IC 


V 


^ C s C m C f 




)( 


C h C H > 

c T c R )J 

(2.5) 


For a safe design, the. constraint on HPW is 
HPW 'y/ HP for all the gear sets. 

The notations used in the above expression had 
already been mentioned in the ’Nomenclature * . 

A brief account of the various factors used in 
arriving at the expression of HPW is given here t 
O verload Factor C p : It evaluates the smoothness, roughness, 
and,*jpeak loads developed by the driving and the driven 
machines. The overload factors used have been taken from 
the tables given in the ASME transactions. 


Dynamic Load Factor C v : It depends upon the accuracy of 
tooth spacing and profile, pitch line velocity, inertia 
and stiffness of rotating elements, transmitted loads, 
tooth and blank stiffness, and properties of the lubricant. 
Commonly used dynamic factors are shown in Figure 2.2. 
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Size Factor C s It reflects the effect of dimensions and 

— — — W 

uniformity of material properties. For normal gears using 
good commercial material; as in this case, C c equals 1. For 
practical design purposes, the necessary adjustments for 
size are usually included with the load distribution factor. 

Load Distribut i on Factor : It measures the effects of 
nonuniform distribution of load across the face width. It 
is the ratio of the maximum to the average load intensity. 
Load distribution depends on - 

1. Cutting errors 

2. Internal Bearing clearance 

3. Run Out 

4. Tooth stiffness 

5. Blank Stiffness and Shaft Stiffness 

6. Bearing Deflections 

The load distribution factor has also been taken 
from ASMS transactions. 

Geometry Factor t Includes the geometric terms developed 
in the formulation of relative radii of curvature, as 
well as the load sharing ratio between the oblique contact 
lines. The AGMA durability standards recommend 

nip 

1 = °« 225 ( 5 ^ ) 

where 

- gear ratio. 



23 


Surface Con dition. Factor CF s Depends on profile surface 
finish, residual stress, and plasticity effects. The 
surface -finish is changed by operation under load, and 
usually stabilises at a definite value regardless of whether 
the initial finish is finer or rougher. Gears generated 
with reasonable care, as well as shaved can be rated with 
C f = 1.0 - 1,20. 

Life Factor C-jU Adjusts for the required number of tooth 
contact cycles. Usually the marine gears are designed for 

Q 

10 cycles/day and the life factor used is 1.10. 

Hardness Ratio Factor : Measures the hardness difference 
between the two meshing gears. Ip the present case, as 
the material of gear and pinion is the same and they have 
been subjected to the same degree of hardness, = 1 is 
recommended for use. 

Factor of Safety Enables the designer either to attain 
high reliability or if desired or necessary to assume a 
calculated risk of failure . Factors which must be consi- 
dered in selecting a safety factor are : Effects of pitting 
on noise generation and uniform motion transmission and 
the possibility of stress concentrations initiating a tooth 
fracture. 

The factor of safety of 1.25 has been used to take 
a. higher reliability against failure. 
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jjl&stic Coeffic ient; CP Depends on the moduli of elasticity 
and Poisson's ratio of the pinion and gear material* It is 
given "by the expression 

Cp = SQRT (1/C71- + ^r~m 

2 .7-5 Bending Strength 

Loads that helical gears can transmit are normally 
limited by the pitting resistance of the tooth profile 
surface. Gear strength (bending considerations) becomes 
critical when materials of high hardness are used and for 
applications subject to high over loads. Though the gears 
will usually fail in wear, it is always advisable to check 
them from their strength point of view* The bending 
strength for a gear set is given as 


-»,<f 787 .040 


hpb = | iBEMi _a_ 

74,000 


r fats 1,26 

v_ n 0.067 


D. 


3 .224 


(^np) 


1.192 


Y 


t* 1.26 


n 


r .2.004 

(Cos^D) 


11 


2.450 


1 


( 2 . 6 ) 


np 


where the symbols used have already been defined in the 
•Nomenclature ' . 

For a 20° pressure angle 


L N 


0.550 - 1.4l/(N n -1.89) 


.756 
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P(f) = f(x) + r y ~~^r 

i=l «!&> 

At any constraint surface i, gj_(X) reaches zero 
thereby shooting up the values of total function p(X) . The 
function P(x) is given as high value as 10 in case any 
one of the constraints is violated. This helps in identi- 
fying the situation that one or more constraints are 
violated. In the algorithm, the solution then starts from 
the previous feasible design. 

2.9 Solution Scheme 

The object is to minimise f(x) subject to N c 
constraints. The constraint list is far from being exhaus- 
tive. Tven with these limited constraints, the solution 
becomes a formidable task because : 

(a) The number of teeth on each gear must be an integer 
i.e., it should be 44 or 45 but not 44.67. 

(b) The availability of cutters restrict the usage of 
standard modules. This is a situation where some of the 
variables can only take predetermined values. 

These conditions make the formulation that of a 
nonlinear problem with integer and real variables. Some 
solutions for pure linear integer programming problems 
have been obtained but the nonlinear part is still 
unexposed. 

The only feasible remedy for this impasse seems 
in floating the variables. After reaching the minima by 
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any of the existing techniques the variables are given 
their next higher integer values. But this method is also 
not free of faults. The final variables must satisfy all 
the constraints. Depending on whether the module value 
is made higher or lower to achieve the next standard value, 
the number of teeth of corresponding gears should be lowered 
or raised respectively. The output speed at the propeller 
must not differ by more than 2% from the required speed. 
Since the rounding-off of the teeth yields a change in the 
transmission ratio, there is a consequent change in the 
output speed. But these may not be as difficult as they 
appear to be. Because of the finite number of design 
variables, the simplest method would be to check every 
possible combination by fixing the range for each variable 
round the mean value found and always checking all the 
constraints. That which leads to minimum area amongst 
the acceptable combinations is the best design. 
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chapter III 

SOLUTION ALGORITHM 

The technique employed in the solution of the 
foregoing constrained minimisation problem is discussed 
in this Chapter. The constrained problem is first conver- 
ted to an unconstrained minimisation problem and is then 
solved as a sequence of unconstrained problems. 

3.1 Conditions for Minimum 

Necessary and sufficient conditions for minimum 
of a function F(X) of n variables is that there exists 
X* for which 

f(X*) f (X) (3.1) 

for all X. If the inequality holds in some neighbourhood 
of X* then the function is said to have a local or relative 
minima at that point. 

The necessary conditions for the relative 
minimum are : 

§4“ = 0 , i = 1,2, n (3.2) 

^ X i 

The set of nonlinear simultaneous equations, thus 
obtained, are not easy to solve. Further, even if the 
above set could be solved there would be no guarantee 
that the obtained solution is not a maximum or saddle 
point rather than a minimum. One way to overcome this 



difficulty is to reduce the constrained minimisation 
problem to an unconstrained minimisation problem. 

3.2 Conversion to Unconstrained Minimisation 

1 constrained optimisation problem can bo cast 
as a mathematical programming problem in the form 
Min f(x) 


g;(X) > 0, 

3 


Subject to ; 

j = 1,2,.....N C (3.3) 

where X is an n-dimensional vector of design variables 
X^, i = l,2,.....n and g^(X) are the given constraints 
on the design. The function f(X) is called the objec- 
tive function and its choice is governed by the nature 
of the problem. 

The constrained minimisation problem is converted 
to an unconstrained minimisation problem as follows : 

L— 


P(X,r) = f(X) + r > 


1=1 


(X) 


(3.4) 


where P(X,r) is the penalty function and r is an arbi- 
trary penalty parameter which in the limit goes to aero. 


The minimisation proceeds over a strictly 
monotanically decreasing sequence of r-values from an 
initial design vector Xq essentially, inside the boundary 
of the constraints and for a positive value of r, say 
(This is what is called 'Interior Penalty Function 
Approach) . The minimum of P(X,r 1 ) must lie inside the 



33 


constraint set because at the boundary some of the g^(X) 
tends to zero and consequently P(X) tends to infinity* 

The minimum of PCX,^) depends upon r , the starting 
value of r, and can 'be written as X(r^)* By reducing 
in the next iteration, the summation term which penalises 
for the closeness to the constraint boundaries is reduced 
and hence in minimising P(X,r) more emphasis is placed 
onF(X). The method has special advantages of not 
following the zig-zag pattern of minimisation (fig*3.l) 
as is observed with other nonlinear programming algorithms. 
This is because the minimisation starts from within and 
does not follow the boundary. 

Fiacco and McCormick^ have shown if 

(1) the interior of the constraint set is non-empty, 

(2) the functions f(X) and g^(X) , i=l,2,....N c are 
twice continuously differentiable, 

(3) the set of points in the constraint sot for which 
f(X) V c is bounded for every finite V c and 

(4) the function f (X) is bounded below for X in the 
constraint sot, then, the optimal solution to the 
unconstrained problem approaches a local minimum of 
the constrained problem as the value of r approaches 
zero. If, in addition, 

(5) f(X) and -g^X) , i = 1,2,....'.N C aro convex functions, 
and, 

P(X,r) is strictly convox in tho interior of the 


( 6 ) 



- constraint set for ever j r> 0, 
then the optimum solution to the unconstrained 
minimisation problem approaches the absolute minimum 
of the constrained problem as r approaches zero* 

3.3 Choice of Initial r 

To start the algorithm the first decision Is to 
be taken about the value of r and the factor C, the factor 
by which the values of r shall be decreased. The algo~ 
rithm, as such, imposes no restrictions on the values of 
r and C except that r ^ 0, C ^ 1. If r is large, the 
function is r asy to minimize, but the minimum may lie 
far from the desired solution to the original constrained 
minimisation problem. On the other hand if r is small, 
the function will be hard to minimize as it then requires 
the initial solution to be very close to the actual 
minima. 


Two methods to' choose the initial value of r 
have been suggested in the literature. 


Cl) r„ 


V f(i 0 )yp(x Q ) 

= 1m> 


(3.5) 


where 






i=l 


gfCvo) 


and represents the gradient of the function 
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Case 1 : ^ 0 

If r^<^ 0, minimisation of I{X.q) alone, without 
considering the penalty term, is carried cut . At every 
new point the r value is checked by equation 3.5. If the 
value Is positive, unconstrained minimisation with 
obtained value of r is carried out. 


Case II t r x = 0 

This means that the unconstrained minimum has been 


reached and 



Case III : r^ 0 

If r > 0, this is taken as the starting value 

for usual minimisation. 

(2) If the initial design is conservative (I.e . , not near 

any constraints) , one would like to pick the initial 

r = r n so that F . (r A ) would not increase drastically 
v min u 

over the original design. In other words, r ought to be 
chosen small enough that in the neighbourhood of the 


initial design, the r terras do not completely 

g ± (X) 

dominate P(X,r 0 ) . A rule which might follow from this 
observation is that if Xq is a conservative design, pick 
r Q so that r Q = is given by 


h = 


fCX 0 ) 



(3,6) 
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This has two weaknesses^ Firstly, the starting 
point may "be too close to some of the? boundaries thereby 
making such selection of r useless, because the r value 
dictated by the above rule might be too small t o allow 
the first minimisation to be carried out. In this case, 
a proper value of Pq will likely be large enough that 
in minimizing P(X,r), F will increase from its value at 
Xq . While this is distressing, it probably cannot be 
helped with this form of penalty function without a. good 
deal of complex logic. Furthermore, unless something 
really happens, very little is lost since r can be 
reduced quite quickly in this method by giving a rea- 
sonable value to C. 

Another approach to this latter problem (initial 
solution being close to constraint surfaces) which seems 
appealing in some cases is to pick a relatively large 
value of r but to temporarily add a new constraint to the 
problem in the form of 

%!0+l = - F «0> } ° 

or to make it easier to -get a starting point 

%c + i ==F « ) - <F( FP + « 5 

where £ is some small amount of increase which will 
theoretically be permitted in F on the first minimiza- 
tion. The penalty function for this revised problem is 
then 

P(X,r)=F(X) + r 
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The minimum for large values "of r is approximately 
the point where the term in brackets is a minimum, ks r 
is decreased, the fictitious constraint term can be 
removed or left in as desired since it will ultimately 
vanish. 


3.4 Extrapolation For r 

The vectors XCr^, X(r 2 ) X (r K ) obtained 

by minimising P(X ,r) over decreasing sequence of 

r- values lie on a trajectory X(r) , where X(Q)is the desired 

solution. From this trajectory, it is possible to 

obtain, b 7 extrapolation, estimates of final solution 

X(0) and the next minimum X( r T< - +2 ) . 

It has been shown the.t the trajectory approximates 
1/2 ( ^ ) 

the polynomial of x . This is because as r tends 

to zero, the function behaves as a linear function of 

1/2 

r . The approx imating function is defined as 

( 2 = 1 ) 

X(r)=a 0 +a i r 1/2 +a 2 r + a K~i r 2 ' ^ 3 * 8 ^ 


where a.’s are the undetermined coefficients. 
1 

For small values of r 

X(r) = X(0) + a*/T 
Solving these for X(0) 

VC" X( r/C) - XCr) 


(3.9) 


(3.10) 


X(0) = _ 

VC -1 

An estimate of minimum point for next r- value 
is obtained bp assuming that 
X(r/C 2 ) = X(0) +aVr/C^ 


(3.11) 
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Solving equations (3.10) and (3.11) for X(r/C^) 
we get, 

X(r/C 2 )= X(r/C) + — (X(r/C) ) - X(t) (3.i2) 

This value of X can be used for the next mini- 
misation, thereby substantially reducing the effort for 
minimisation of POO . However these extrapolated values 
have to be first checked so that they do not violate the 
constraints. 

Estimates of X(0) become reasonably good after 
3 or 4 trials, and towards the end offer more accurate 
estimation of values of X(0) . 

3.5 Method of Unconstrained Minimisation Davidon- 
Fletcher- Powell Variable Metric Method 

This method completely avoids the need for 
evaluating second derivatives and performing matrix 
inversions and yet the sequence of itera.tes converges 
quadrat icaliy to the minimum point, X m . 

Tho method is based on the properties of a 
quadratic function and is designed so that when applied, to 
a quadratic, it minimises the function of n variables in 
’n’ iterations. 

Central to the method is a symmetric positive 
definite matrix which is updated at each iteration. 

It collects in itself the information about the 
curvature of the function P(X) and supplies the current 
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direction of move dj_ by multiplying it with the current 
gradient vector G^. 

A typical iteration in the minimisation of 
F(X) proceeds as follows s 

(1) Given the starting point X^ and the gradient of 
P(X) at X^, dj_(X), the direction of (i+l) t,a minimisation 
of P(X) is given by 

\ • G i (3.13) 

(2) Find <r\ so that FCX^^ df) is minimum along the 
line dj_ . 

(3) Set X 1+1 = X x + cTi* h- 

(4) Calculate the new gradient vector 


Vl = 

and set 



(5) Calculate the new ^H] matrix, Pi + iP y 

, a. a, 1 rH,1 y, Y- "I'h =7 

t 'u 1 = fH 1 + * _ A , A ^ isL-AJ i — 

H i+iJ l H ij + ri Ti - ~ 


h *± 


V&3 h 


(3.14) 


(6) Begin the next iteration from (l) 

The basic algorithm is extremely pove rf ul for a 
first order method, converging quadrat ic aL ly and posse- 
ssing very good stability. By stability j moan that even 
in highly distarted and eccentric functions, it continues 
to progress and needs little of the sort of special 
attention as required by its parallel algorithm, the 
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conjugate gradient method. There is a plausible argument ' 

for this increase in stability in that with conjugate 

gradients, the entire history of the path is carried to 

"kks intelligence of £. d. , a single vector 

2 / a * -3- X 

(j$i ~ | | ) • In this variable metric method, 

Wo carry the data in a full matrix which we carefully 
upgrade at each iteration. Another advantage of this 
method is that in conjugate gradients the carry over term 
^3-d^ is only good If applied to and produces 

non-sense if applied to gradient at some other point, 
on the other hand it can be shorn that £h.Q, is a positive 
definite approximation to the matrix of second order 
partial derivatives and is applicable in the whole space. 

There are certain drawbacks in the algorithm 
which must be kept in mind while applying this method. 

The positive definiteness ofQiJmatrix is preserved in 

)jc 

theory only If <f. is the true minimum point (I.e. 

T 1 

G. _ _ d. = 0) and furthermore, round off error will 

1+1 * i 7 

again dog our steps so that the process can get into 

trouble • 


In applying the method, therefore, care must be 

exercised to ensure that the (ff3 matrix is not updated 

* 

with data arising from poor approximations to (f\ . 

There are a number of approaches to this problem : First 

the algorithm used for computing <jn may be roapplied 
T 

until d. * G* is sufficiently small, another alternative 
is simply to skip the update cycle £_step 3 in the 
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algorithm} when is too large. In other words if 

<fi ls not close enough to minimum along dj_, set H i+1 =H 1 

and _ d i + l ~ “ H ’i+i* G i+i ^- nd continue as before. As long as 
^ ^i +1 ) ^ j the method will continue to progress, 

towards the minimum. 


It is difficult to choose between these approaches, 

the first may require excessive computation to refine 

at points far from X m while on the other hand, the second 

approach may pass up valuable opportunities to improve the 

pi}matrix. A reasonable compromise is to set a moderate 

T 

criterion for dj_ * limit the number of refits for 

<0FI and then skip the update cycle If the criterion is 
not met on the second try (This logic has been used in the 
programme) . 

3.6 The Gradient 

At every iteration the gradient has to be 
calculated in the minimisation algorithm. It can best be 
determined by writing down the exact derivative of PCX) 
with respect to the design variables. Because of the 
large amount of work involved In calculating the required 
partial derivatives, it has .been decided not to derive 
analytic expressions for the components of the gradient 
vector, but instead, to approximate the gradient by 
difference equations which involve only function evalua- 
tions. 

The gradient of the function P(X) with respect 
to the design variables Is given by* 
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G(X) = 7«i) = (^1 ^JE., ~ap vi 

It may also be written as 

N 

ps, = ( Hl- r gij. 

i=l 1 1 


(3,15) 



(3.16) 


The variable metric method technique for minimisation is 
based upon the assumption that the function being minimised 
can be satisfactorily represented by a quadratic function 
in the vicinity of the function minimum. Therefore 


expanding the function F(X) of one variable about point 
X 0 by Taylor series up to quadratic terms we have 


“Xq )"*■*«# 

(3.17) 

Neglecting the second order and higher order terms 
dF* _ F(X X ) - F(X 0 ) 

dxt^ ar-~V” (3.18) 

Close to the minimum, the function P(X) has very 
large curvature indicating significant values for the 
second and third derivatives which have been neglected 
in the forward difference scheme. Even if g(x) and f(X) 
are nonlinear, still they are quite well behaved as 
compared to P(X) at any point inside the domain. Thus the 
gradient from the equation (3.16) should be used. 


P(X X ) = F(X 0 ) - (g) 


X 


0 


(x i • V + 1 
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3.7 Linear Minimisation (Fibonacci Search) 

The one -dimensional search procedure used in the 
programme is the method called the Fibonacci Search. Here 
we make use of the Fibonacci Numbers, which are defined by 

%*•! = f N + r U-l > F o = h = 1 
The principle of the working of this method is as follows : 

Assume that the minimum of F(x) along x lies between 
two points A and B, According to the accuracy of the 
result required, choose a value of N , such that 

X A ~ X B , r 

f n ^ x 

Choose two points 1 and 2 at distance scK^ and rX. 2 f r0Iil A. 

oCi = *c a = 

Compute the function values at 1 and 2 as F^ and Fg. Assuming 
that the function is unimodal, i.e. the point with smaller 
function value to be nearer the minimum compare F^ and Fg. 

The portion outside the larger value is discarded, and the 
search is repeated within the rest of the zone, until the 
minimum is reached to the given accuracy (Fig. 3. 2). 


3.8 Slope in d. Direction: 

In one dimensional minimisation of P(X) , slope 
in the direction d is required. This can be obtained in 
two ways : 
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1. By taking dot product of the vector d and G at the point 

\f 1 = d T • G 

2. By taking finite difference in d direction. 

the fiist method, to calculate the gradient, it 
is essential to e valuate the function value n times. 
Whereas, in the second method only one computation of the 
function is necessary for the forward difference scheme. 
Therefore the second method has been preferred here. 

The forward difference scheme in the d direction 
is as follows i 


= *± + e 


a i 


where f- is a small step in d direction. The slope tyKcT ) 
in d direction is then 

d V (<r ) _ df(e ) 

de 




t _ 


ae 


- r 


N c 

V dgi 1. 

d£ e 2 
i=l 6 i 


(3.23) 


whe re 



df(X) = f ( 
d<? 

dgj s 
d£ 

3.9 Convergence Criteria 

The following are the various convergence criteria 
used in the algorithm at various stages* 

1* For one dimensional Minimisation 


a) 


— T — 

d - G 


d * G 


< e si 


(3.24) 
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where £ is the cosine of angle between d and G. 

b) The percent change in the function is less than specified 
tolerance 

4-£s2 (3 - 2S) 

2. For minimisation for a particular r 
a) When the approximation for the percent change in P(X) 
in the next iteration has become less than a certain 
pre-assigned -value. 


f(X) 


GUIS ^ (3,26) 
PCX) ^ ml 

b) Another criterion to stop the iterations for a parti- 
cular r is when n/2 consecutive one -dimensional minimisa- 
tions have yielded (step in the d direction) less than 
one percent of the smallest value of the design variable. 


3. For the Unconstrained Minimum 

a) When the percent change in the objective function 
between minima corresponding to two consecutive r values 
is less than the given tolerance, that is 


-min 


(X, r-l) - f. 


min 


(X,r) 


f min <r) 


<0 £ 


u2 


(3,27) 


b) When the minimum percent change in any design variable 
corresponding to two consecutive minimum r values is less 
than the given tolerance. 
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X^r/C) - X . ( r) 

~ r T 



m 


C) When the percent change in the extrapolated 
value f(0) is less than the given tolerance. 


( 3 . 28 ) 

function 




M|§Sfl| pipi 

■sift 

N|0IM$ 

MSM 

■MMi 

B ill Hw 

HMlilf | :: ' * 

1 


HBHHBi 

MH8I 


■■ 









49 


CHAPTER IV 

RESULTS AND DISCUSSIONS 

Three illustrative examples have been solved to 
show the applicability of the unconstrained minimisation 
algorithm to solve the marine gearing problem . First 
example is the main reduction gearing of a- Naval 
destroyer. Results obtained therefrom are tabulated. In 
the second example the gear box for a cargo ship - JAL 
DEEP, having 60,000 K.P. has been designed and the 
results compared with the existing design* In the third 
example the gears have been designed to transmit 75,000 
H.P., but the results could not be compared as the exis- 
ting results for this were not available* 

In all the previous three examples, the speed 
reduction had been from 3000.0 R.-P.M. to 150.0 R.-P.M., 
the gear ratio being 20.0. The criterion for ‘best 1 
design in all the discussed examples is the minimum deck 
area occupied. Constants which are common to all the 
three examples have been tabulated in Table 4.1. Table 
4.2 shows different initial starting vectors for the 
discussed example s.- 

4.1 Example 1 

This is the main drive for a Naval destroyer 
having 40,000 H.P. at its propeller. The drive has two 
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reduction stages, each stage employing double helical 
gearing. Input information for this problem is given 
in Table 4.1 and 4.2. Table 4.3 compares the existing 
and optimal results. Table 4.6 shows various optimum 
solutions. 

4.2 Sample 2 

This is the reduction gearing for JAL DEEP, 
a cargo ship with 60,000 H.P. at its single screw. The 
speed reduction is carried out in two stages so as to 
facilitate the accommodation of bull wheel in the space 
between the centre of the shaft and the bottom of the 
ship. Double helical gears have been used. Table 4.4 
compares the existing and the optimum results* Table 
4.7 shows the variation in the optimum solution with the 
variation of the starting value of ’R f and different 
starting vectors. 

4.3 Example 3 

This is for the drive of a single screw ship 
with 75,000 H.P. The speed reduction had been carried out 
in two stages using herringbone gears. Table 4.8 shows 
the final optimum results obtained by the variations 
of 'H 1 and starting vectors. 

4.4 Discussions 

The results ' .from the first example show about 
23.4% reduction in area from the existing design. The 
results of second example also show a reduction of 54.3% 
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in area from the existing design. The rounded numbers of 
teeth have adjusted themselves in such a fashion that out 
of two mating gears, the numbers of their teeth have no 
simple common factor, which is desirable from the wear 
considerations. This added advantage shall be obtained 
almost always because the transmission ratios in both 


the reduction stages do not arrange themselves in simple 
fractions. The same objective could also be achieved by 
adding a hunting tooth on the gear, in case there is any 
common factor between the number of teeth on gear and 
pinion. 

The objective function in the first example 
decreases only by 23.4$ whereas in the second example, 
the reduction is as high as 54.3$ . This is so because 
the first example was for a Naval ship wherein the utmost 
attention had been paid to save as much spaco as possible, 
whereas, the second example is for a cargo ship where 


different criterion, other than space , must have 
dominated in the design philosphy. The optimum results 
for example 3 could not bo compared duo to unavailability 


of the existing results. 

The percentage improvement in the objective 
function in the discussed examples, may be attributed to 
one or more of the following points. 

1. The method of existing gear design may be different 


from the one used here. 


i i. I.T KANPUR, 

I CE.-; T AMT, 




2- In the present algorithm, the computer facility has 

been fully exploited to find out the best combination of 

the independent design parameters to achieve the optimum 

design, m conventional design procedures, this is not 
possible . 

3. Better material with 5S0 BHN has been used in the 

present design compared to conventional gear materials 
with BHN of 400-500. 

The optimum search is possible here because the 

algorithm seeks the minima fat least local) in the feasible 
design space. 

The author docs not claim that the global minima 
had been reached in the examples discussed. This is 
because firstly not many starting vectors were chosen 
(only 3 starting vectors for each example) due to shortage 
ox computer time and even for one starting point the 
iterations were stopped when no significant changes 
occured in the objective function. Secondly even in 
some of the examples with a particular starting vector, 
the convergence was rather slow a s compared to other 
starting vectors for the same example. This may be due 
to the inefficiency of the linear minimisation algorithm* 
which prevented the process from making satisfactory 
progress towards the optimum solution. It Is found that 
when the starting point of the linear minimisation lies 
close to a constraint boundary, so that even a very small 
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step away from this point causes a constraint to be 
violated, the linear minimisation procedure comes up 

with a poor approximation to the minimum of the function 
along that line. 

Another problem which was encountered in using 
the present method is the choice of an initial value of 
penalty parameter *R* . m Chapter 3, two methods of 
selecting the initial value of *R« have been discussed, 
both of which depend upon some knowledge of the charac- 
teristics of the gradient vector. In using the finite 
difference formula to compute the gradient at each step, 
these necessary facts are not available. It would be 
better if a more rigorous method for choosing the startin, 
value of »r« could be found. 
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TABLE 4.1 

VALUES OF CONSTANTS COMMON TO ALL EXAMPLES 


SYMBOL IDENTIFIER 


O' 

PHIA 
RPM l 
RPM 3 

E 

rr 

EG 

6- 

AK 

C 3 

C3 

Ug 

UG 

Up 

UP 

SALE 
FATS 
SHEAR 
ITIR 
, TRAN 
EPS1 
EPS 2 
EPS3 
EPS4 

c 0 

CO 

c r 

CR 

Cl 

CL 

C s 

CS 

Ch 

CH 

Ct 

CT 

Cf 

CF 

SENS 

PMOD 


MAGNITUDE 

20 . 0 ° 

3000.0 

150.0 

25000000 .0 lb/in 2 

861.00 
1.67 
0.20 
0.20 

225.000. 0, lb s/in 2 

150.000. 0, lbs /In 2 

850.000 lbs/in 2 
200 

20.0 
0.0005 
0.0005 
0.0005 
0.0005 
1.50 
1.25 
1.10 

1.00 

1.00 

1.00 

1.00 

0.50 

1.0, 1.125, 1.375, 1.50, 

1.75, 2.0, 2.25, 2.50, 

2.75, 3.00, 3.25, 3.50, 

3.75, 4.00, 4.5, 5.0, 

5.0, 6.0, 6.5, (.0, 

7.5, 8.0, 8.5, 8.0, 

9.5, 10. 0, 11.0, 12.0, 

13.0, 14.0, 15.0, 16.0 


Column of symbol is kept clank where the 
symbol and the Identifier are the same. 



different initial di 


Representation X = 


> x- 


EXAMPLE 1 : 

50.0 
80.0' 

2.40 
.2.40 

25.0 

X = <J25.0 > 

1 7 5.50 i 

55.0 

8.0 

9.50 

. 0.00065 
EXAMPLE 2: 

* 100. 0 

130.0 

2.0 

2.0 

, 28.0 
< 28.0 
^ 5.0 

70.0 

8.0 

8.0 

0.00085 
EXAMPLE 3 : 

150.0 

150.0 

2.0 

2.0 

30.0 
Xn= 30.0 

■ 5.0 

65.0 

8.0 

8.0 

10.00085 J 


N.B. - The subscripts 
and 3 . 


TABLE 4.2 

iSIGN VECTORS FOR EXAMPLES 1, 

'ZPl 
ZF2 
PNDl 
PND 2 
SIA1 

< SIA2 > 
i TRANl < 

ALPHA 
GNFI 
GNF2 
ERR 


bo 


2 AND 3 


X 2 = < 


' 45.0 

80.0 
2.45 
2.45 
24. C 

24.0 

5.50 

52.00 

8.00 

8.50 
0.00070 

75.0 

100.0 
2.20 
2.20 

26.00 

26.00 

5.00 

65.00 

8.00 

8.00 

0 .00080 


> 


> 


> n = 


'41.62 
86.65 
2 .314 
2.312 

„ , 22.81 
X 3 = <T 22.86 
5.87 
52.27 
10.90 
11.55 
0.00055 


42.0 

86.0 
2.20 
2 ,3i5 
22.80 
22.80 
5 .89 
52.50 
10.00 
11.55 
0.00055 


Xs =< 


f 

> 


41.62 


86.60 


2.314 


2.312 


22.81 


22.86 

5.87 

? H-t 

52.27 


10.90 


11.56 

!■ 

0.00055 . 



1,2 


60.0 
120.0 

2.50 
2.45 
25.00 

25.00 

6.00 

55.00 

10.50 

11.00 
0.00055 

1,2 5 3 stand, for starting vector 
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TABLE 4.3 

COMPARISON OF OPTIMUM DESIGN TO THE EXISTING 
DESIGN FOR EXAMPLE 1 


3. No 

. K description 

$ 

^EXISTING DESIGN . ^OPTIMUM 
5 1st $ ^2nd 5 1st ~ 

k Stage 1! Stage 5 Stage 

DESIGN 

5 2nd 

8 State 

1 . 

Number of teeth on pinion 

58 

115 44 

80 


Number of teeth on gear 

2 SO 

460 259 

281 

2 « 

Module 

7.8* 

7.8* 10.0 

10.0 

3 • 

PCD of pinion in m.m. 

452.4 

897.0 440.0 

800.0 


PCD of gear in m.m. 

2262.0 

3588.0 2590.0 

2810,0 

4. 

Helix angle in degrees 

27.93° 

27.93° 22.90 

22,90 

5 . 

Transmission ratio 

o.O 

4.0 5.78 

3,46 

Obtained by converting the circular pitch 

into equivalent 

module . 

N.B. 

- The Optimum Design refers 

to Initial Design Vector-: 

2 


TABLE 

4*4 




COMPARISON OF OPTIMUM 

DESIGN TO 

THE EXISTING 



DESIGN FOR EXAMPLE-2 




1 ' " " TEXISTING DESIGN ii OPTIMIST DESIGN." ' 


S .No 

.5 DESCRIPTION 

5 

5 1st 
| Stage 

ft 2nd 

8 State 

5 1st 

5 Stage 

5 2nd 

5 Stage 

1. 

Number of teeth on pinion 

46 

97 

34 , 

76 


Number of teeth on gear 

255 

354 

207 

260 

2. 

Module 

11.82* 

10.74* 

13.00 

3.1.00 

3. 

PCD of pinion in m.m. 

543.80 

1042 .0 

442.0 

836.0 


PCD of gear in m.m. 

3020.0 

3813.0 

2691.0 

2860.0 

4. 

Helix angle in degrees 

28.75® 

28.75° 

23.08° 

23.08° 

5. 

Transmission ratio 

5*55 

3.64 

5.94 

3.37 


* Obtained by converting the circular pitch Into equivalent 
module . 

N.B. - The optimum design refers to initial design vector-3 



TABLE 4.5 
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OPTIMUM DESIGN 

FOR EXAMPLE-3 



S « 

T~ 

8 DESCRIPT- ION 

£ 

it JOFTIMUM DESIGN 

l EXISTING DESIGNS 1st T 2nd 

1 . 5 Stage 5 Stage 

i. 

Number of teeth on pinion 

Results not 

' 43 

87 


Number of teeth on gear 

available 

255 

302 

2. 

Module 


11.00 

11.00 

3 * 

PCD of pinion in m.m.. 


473.0 

957.0 


PCD of gear in m.m. 


2775.0 

3322.0 

4 • 

Helix angle in degrees 


22.93° 

22.93° 

5. 

Transmission ratio 


5.82 

3.44 

N.B. 

- The optimum design refers 

i to initial desi 

gn vector 

-1 


TABLE 4.6 


OPTIMAL SOLUTIONS WITH DIFFERENT INITIAL DESIGN SECTORS 

FOR EXAMPLE -l 


Star- $ 

FIRST OPTIMAL 

Tsecond 

OPTIMAL 5 THIRD 

OPTIMAL 

ting 5 

SOLUTION 

8 SOLUTION 

l SOLUTION 

V alue l 

Area in 

5 Time in 

l Area in 

igqmgKHi 

5 .Area in 1 Time in 

of R 8 S a .Meters 

5 Min-Se c s * 5 S a •Me t e r s in- Se c s . 

xSqJMete 

rs$Min~secs. 

FO/FP, 

1.48 

11-34 

1.50 

7-16 

1.47 

6-01 

10.0 

1.44 

7-51 

1.40 

6-46 

1.44 

7-50 

50.0 

1.46 

6-53 

1 *43 

6-55 

1.48 

6-10 

100.0 

1.51 

6-43 

1.44 

6-19 

1.46 

6-47 


Value of optimum solution = 1.40 sq.meter 
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TABLE 4.7 


OPTIMAL SOLUTIONS WITH DIFFERENT INITIAL DESIGN 'VECTORS 

FOR EXAMPLE- 2 


Star- 5 FIRST OPTIMAL $ SECOND OPTIMAL ~‘l “ THIRD tipTIMAL ” 

ting L_. SOLUTION l SOLUTION l SO LUTION 

Value 3 Area in T Time in \ Are a~£rTT T ime in f Are a ~in T T~ime“ in 
of.., £ 1 S a_.Mete i^iMin-se c s . ESa.Meters jlHin-secsiS a .Meters Illih-secs . 


FO/FP 

1.69 

10. 0 

1.66 

50.0 

1 .66 

100. 0 

1.55 


7-10 ' 1.65 

6-41 1.58 

6- 44 1.5S 

7- 11 1.59 


6-14 

1.58 

5-36 

6-12 

1.50 

6-21 

7-08 

1.72 

7-04 

7-13 

1.54 

6-09 


Value of optimum solution = 1.58 sq.meter 


TABLE 4.8 


OPTIMAL SOLUTIONS WITH DIFFERENT INITIAL DESIGN VECTORS 

FOR EXAMPLE-3 


Star- 5 FIRST OPTIMAL 
tins & SOLUTION 

T SECOND 

5 SOLUT 

OPTIMAL 

'IONS 

T THIRD OPTIMAL 

5 SOLUTION 

Value 5 Are a in 

8 Time in 

Q Area in 

5 Time in 

fi Area in 

£ Time in 

of R jjSq.Meters 

| Min- secs 

. jjSq. Meters 

; jjHin-secs 

. jjSq. Meter 

|M in- secs. 

FO/FP 1.87 

4-55 

2.03 

5-47 , 

1.98 . . 

.. 7-28 

10.0 1.96 

6-44 

2.11 

■ 6-41 

2.08 

6-46 

50.0 2.03 

7-07 

2.09 

7-01 

2.03 

7-68 

100. 0 1.95 

7-23 

1.97 

6-25 

2.10 

7-21 


Value of optimum solution = 1.87, s.q,.me ter.. 



comparison of results 


EXAMPLE j : 

Value of the objective function 
of the existing design 

Value of the objective function 
of the optimum design : 

Percentage Reduction in the 
Deck area ' 

EXAMPLE 2 : 

Value of the objective function 
of the existing design = 

Value of the objective function 
of the optimum design = 

Percentage Reduction in the 
Deck area 

EXAMPLE 3 s 

In the absence of the data for this 
the results could not be compared. 


: 7*6 sq.m. 

5,82 sq.m 

7.6 - 5.82 
7.6 

23.4$ 

10.30 sq.m. 

5.70 sq.m 

10 .30-5.70 

10.30 

54.3$ 

example , 
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CHAPTER V 


C owe IX 


Slows AND ESC 0 MMEND AT IONS 


The application of the proposed algorithm to 
solve a marine gearing^ problem shows that the c duplicated 
machine design problem can be tackled by computer.- Whether 
the algorithm is of any practical use or not depends 
upon the problem and the situation in which it is 


encountered. In general, the question as to whether a 
mathematical programming approach is suitable or not 
depends upon the following considerations t 

1. The need for obtaining an optimum design.- 

2. The effort and cost required to write a programme 
identifying the variables and constraints in the 
problem compared to manhour requirements in a conven- 
tional design office for the same problem, 

3. The capability of the programme to tackle a wide 
range of similar problems. 

The programme , that lias been developed, is very 


general and can handle any typo of marine gearing design 
problem to suit the various requirements of the designer. 
Therefore it spreads out the cost of developing the 
programme and the computer time. As a result, the 
economy due to the computer aided design procedure in 
terms of time and engineer-hours in comparison with the 
conventional design procedure may be more* 
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The author does not claim to give sufficient 
information for final manufacturing of the gears, but 
then, the results of the programme definitely allow the 
designer to draw the final drawing and obtain other 
important dimensions from the assembly drawing itself* 

Thus the routine part of the analysis is not. needed and . 
using computer the designer obtains an optimum design 
as a bargain. Finally, it can be stated that the time 
and cost of developing the general programme, though 
considerable is quite insignificant to those for preli- 
minary design of even a single stage reduction gearing 
to tailor the various geometrical and strength require- 
ments. The added advantage of programming is that a 
number of optimal designs with different degrees of relia- • 
bility can be obtained with least efforts. This makes 
room for the designer to use his discretion to pick up a 
suitable design. 

5.1 Recommendations 

The author recommends the following as a basis 
for further work , in the field* 

1. One very severe short coming of the algorithm is 
that it is not self starting. One has to feed starting- 
feasible solution. However, remedying this is not as 
difficult as it sounds. The initial number of teeth and 
module can be given very high values thereby leaving no 
room for inadequate power transmission. But the 



difficulty is likely to come in satisfying the various 
geometrical constraints. With few hit and trial efforts, 
one can find a feasible solution. 

2. The programme could be made more general by inc orp o rating 
the design of gear box. for multiscrew ships. However, this 
requires an entirely different layout of the reduction 
gearing, but the same could be achieved from the present 
programme without much efforts compared to the amount of 
versatility the programme is likely to acquire. 

3. One of the difficulties encountered was the erratic 
behaviour of linear minimisation algorithm at points near 
the constraint boundaries. It is suggested that a better 
way of improving the method may be to modify the linear 
minimisation routine, and tailor it specifically for 

the particular example being considered rather than 
making it more general for all types of problems. 

4. In gradient calculations, the derivatives of cons- 
traints which do not vary with respect to certain design 
parameters have also tc be computed for convenience in 
programming, though it is obvious that the derivatives 
would be zero. This accounts for much computer time. 

An. algorithm may be developed in a minimisation programme 
which would identify those constraints which do not vary 
with respect to a variable about which the derivative is 
being calculated and skips the derivative computation. 
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CORRECTION OF GEAR TOOTH 

The number of gear teeth which have been obtained 

from the optimisation programme shall have- non- integer 

values in general* They must be rounded off to the next 

integer with a view that the final speeds at spindle do not 

change by more than 2 %, From strength considerations, if 

the pinions are given positive corrections, their strength 

increases at the cost of already overdesigned gears. Total 

correction required to accommodate the modified number 

of teeth is .given as 

( INV( oQ -INV( . 'q ) ) « (Za* + Zb 1 ) 

X T 2 tan (o(o) 

where Za' = modified number of teeth in pinion 
Zt) 1 = modified number of teeth in gear 
= modified pressure angle 

= tan 1 
Z 1 + Z b T 

and Q = Za ^Z “ * tan C< 0 

where Z = unmodified number of teeth in pinion 

a 

Z b = unmofied number of teeth in gear. 
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APPENDIX B 
COMPUTER PROGRAMME 

This Appendix gives the descriptions of nine 
subroutines which add up to make the optimised design of 
marine gearing. A listing of each subroutine is included* 
The input and output formats are described and a part of 
output listing is printed. 

B.l The Programme 

The programming has been split in. two phases, 
optimisation and design. The first part is carried out 
in the main routine which requires another six subroutines 
to reach its end. Second phase is carried out in the 
subroutines Modif 5 Design and Corree, which needs the 
optimum design vector from main as input information. 

The purpose of different subroutines together 
with the necessary input for them follows in the subsequent 
paragraphs. 

B.l-l Main ; 

Purpose - To read the input data to extrapolate the 
design vector for a new starting 
point (to be used in subroutine 
Uncons.) and to stop iterations when 
convergence criteria are satisfied. 
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B.l-2 Subroutine Const. : 

Purpose - To calculate the various expressions used 
to calculate the constraints in the sub- 
routine 'Penal* , and to design the gear 
shafts* 

B.l-3 Subroutine Pen a l(F0 ,FP,PT ,X,N,R,GX) : 

Purpose - To calculate the value of the objective 
function and to evaluate the nomalised 
constraints* To set the value of penalty- 
function as 10 3 ° whenever a constraint is 
violated . (This identifies an acceptable 
solution from the unacceptable one). 

Input - The present design vector, number of variables 
and the value of r in the unconstrained 
minimisation are the inputs to this sub- 
routine. For every evaluation of constraint 
the subroutine first calls the subroutine 
’Const* to modify different expressions, 
used in the constraint evaluation. 

B.l-4 Subroutine UnCons. (X,H,R.) : 

Purpose - To carryout unconstrained minimisation for 
any ’ r' value in the sequence of r mini- 
misations. In process, to identify a 
situation as whether to modify H matrix 
or not and to stop iterations when the 
convergence criterion are satisfied. 
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Input - This subroutine is called from Main every 
time when ’r* is reduced. Starting design 
vector, number of design variables, and the 
value of r are the input information. This 
subroutine uses three other subrout in es 
namely PENAL , SEARCH, HMATRI at different 
stages. Convergence limits are transferred 
from Main through common statements. 

B.l-5 Subroutine Search (STEP, X,R,N) 

Purpose - To carryout the one dimensional m in imi sat ion 
and evaluate a step size which makes the 
penalty function minimum in the present 
direction of move. This subroutine is called 
from Main and ^nconst 1 each time a new 
direction is found out. 

Input - The new direction and the present v al ue of 
design vector are input here. Convergence 
criterion has been permanently defined in 
the subroutine itself. 

B.l-6 Subroutine Grad (G,X,N,l) ; 

Purpose - To calculate the gradient at any point in 
the design space by forward interpolation. 
This is called from TJncons every time a 
step in the old direction Is taken and a 
new design vector is reached. 
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input - Inputs are the present value of X and number 
of design variables. The step length used in 
the forward difference scheme is permanently 
defined in the subroutine, 

B.l-7 Subroutine HMATRI (H, SIG, h) : 

Purpose - To modify the Hessian matrix used to calculate 
the direction of move at every point. Every 
time step size in the previous direction 
meets the convergence requirements, this 
subroutine is called from Uncons. 

Input — Inputs for the subroutine are the present 

Hessian matrix, the value of the step in the 
previous move, the difference of the present 
and previous gradient vector and the number 
of variables. 

B.l-8 Subroutine MOD IE (X,N, R) : 

Purpose - To round off the number of teeth of different 
gears and pinions to make them integers. To 
round off the module of gears to the next 
standard module available. 

Input - The final optimum design vector. 

B .1-9 Subroutine CORREC (ZP, ZG, TRNS, XT, XP, XG) : 

Purpose - To distribute- the total correction on the 
pinion and gear for achieving the compati- 
bility. The information about the graphs in 
I. S. (3756) used for this purpose is permanently 
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defined in the subroutine in data statement. 
This subroutine is called for correction of 
each pair of gears. 

Input - The total correction value obtained from the 
subroutine DESIGN, rounded number of tee th- 
in gears and pinions and the transmission 
ratio between them. 

B.l-lO Subroutine DESIGN (X0 S N, NC) s 

Purpose - To round off the centre distance between 

the gears. To calculate the total correction 
for different sets of gears and to convert 
EPS system into MKS system. 

Input - The final design vector from subroutine 
MOD IF . 

B .2 Programme Output : 

The input data is written out in the same format 
in which it is read. 

Other outputs are presented with self explaining 

headings. 

B .3 Programme Listing and Sample Output s 

The subsequent pages contain the listing of all 
the subroutines and a portion of output. 
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16 

23 

25 

27 

28 
29 

31 

32 

33 


34 


35 

36 

37 


MAIN 

DIMENSION XMIN(15), XN(15), SPAN(15), Xl(l5\ 
DIMENSION GX(4Q), GXl(4o), GX2(4Q), D{15), 0(15) 
DIMENSION X (.1 5 ) » XR ( 15 ) , XL ( 1 5 ) , XO(15), G(15) 
COMMON/ DMSl/ PI, RAD, NC 

COMMON/ DMS2/ COl, CR1, CL1, CS1, CHI , CT1, CF1 

C0MM0N/DMS3/C02,CR2,CL2,CS2,CH2,C72 ,CF2 

COMMON/ DMS4/ HP, PHI AN, EG, UG » RF Ml * RPM2, RPM3 

COMMON/DMS5/DPl,DP2,DG.l,DG2,31,B2,ZGl,ZG2,PLVl,PLV2. 

COMMON/ DMS 6/ Cl, C3, AK, GO, CF, SENS, N 
COMMON/ DMS7 / SALBl, SAL32, FATS, SHEAR 


COSB 


COMMON/ DMS 8/ CDl , 
COMMON/ DMS9/ TORI, 
COMMON/ DMSl 0/ EPSl, 
COMMON/ 

COMMON/ 

COMMON/ 


0ZG2 


SC2 


CD2, CDS » QZP 1 » 0Z61 , 0ZP2, 

T0R2, T0R3, bMI, BM2 , P,M3 
FPS2 , EPS3, EPS4, ITIR 
DMS 11/ ZPl, ZP2 , PND1 , P N D 2 , SIAD1, S I A D 2 
DMS 12/ SlAl, SI A.2 , TRAN1, TRAN2 , ALPHD, ALPHA 
DMS 13/ CNF1 , CNF2, ERR, DSH1, DSH2 , DSH3 , SCI 
COMMON/ DMS 14/ XMIN, XN, SOAN 
COMMON/ DMS 15/ GX , GXl, GX2 

COMMON/ DMS 16/ HPW1 , HPW2, +PE1, HPB2 , HEIGHT, EASE, AREA 
C0MM0N/DMS18/ CORRPl, CORRG n , C0RRP2, CORRG 2 
COMMON/ DMS 17/ D, Y, G 

COMMON/ DMS 20 / FCR1, FCR2, PCR1, PCR2, WWl, WW2 , WDi, WD2 
COMMON/DMS 19/ PMODl, PM0D2, TRAN, FOB 

FORMAT! /, 10X, * EXTRAPOLATED VALUES OF VARIABLES AT NEXT 
1 R- VALUES*,/, 3 ( 10X, 5pi0o6 / )) 

FORMAT ( ///, 1 OX , * NORMALISED OPTIMUM DESIGN VARIABLES* , / , 

1 3 ( 10X, 5Fl0 o 6 / )) 

FORMAT( /, 10X, * ACTUAL' DESIGN VAflABLES*, 

FORMAT ( / 10X»* DIAMETER OF FIRST PINION®*, 

1R OF FIRST GEAR = *,Fl2„6 ) 

FORMAT ( / 10X>* DIAMETER OF SECD, PINION®*, F 12 0 6 , / / 1 OX ,* DIAME23 
1R OF SECD 0 GEAR=*, 612,6 ) 

FORMAT ( / 10X ,* THE FACE WIDTH OF FIRST GEAR SET=*»F12,6, 

1 //10X»* THE FACE WIDTH OF SECDoGEAR SET = * , FI 2 • 6 ) 

IOX,* THE PITCH CONTACT RATIO OF FIRST GEAR SET=*»Fl 2 c 6 

Ra.tio of first gear set=*,fi2 0 6 > 

CONTACT RATIO OF SECDo GEAR SET**,Fl2 5 S 
RATIO OF SECDo GEAR SET=* ,F12 .6 ' ) 

LINE VELOCITY OF FIRST GEAR SET=*,F15.6 
1» // 10X,* THE PITCH LINE VELOCITY OF SECD 0 GE*-R 


/, 3 ( 10X,5Fl0o4 / ) ) ) 
F12 e 6 ,//10X»* DIAMETE 


FORMAT ( / 

1,// 10X»* the 
FORMAT ( / 10) 
1»// IOX,* THE 
FORMAT ( / 10X, 


9 ' 

FACE CONTACT 
:»# THE PITCH 
FACE CONTACT 
THE PITCH 

// iox. 


•V. 


2 SET® # » F 1 5 o 4 ) 

FORMAT ( / 10X»* T8E MAXIMUM HoP OF THE FIRST GEAR SET ON WEAR Bl 
1 SiS=*,F15o 4,// IOX,* THE MAXIMUM HoFo OF THE FIRST GEAR SET ON BtlH 
2D ING STRENGTH BASIS®*, F15 b 4 j 

FORMAT ( / IOX,* THE MAXIMUM HeP OF THE SECDo GEAR SET ON WEAR 
1SIS=*,F15 0 4,// 10X»* THE MAXIMUM HoF e OF THE SECDo GEAR SET ON B4^ 
2DING STRENGTH BASIS®*, Fl5„4 j 

FORMAT ( / IOX,* NUMBER OF TEETH OF FIRST P I N 1 0N=* , Fl2 0 6 ,// lox , *N 
1UMBER OF TEETH OF SECDo PINION®* , Fl2 0 6 > 

FORMAT ( / IOX,* NUMBER OF TEETH OF FIRST GEAR =*,F12 0 6,//lOX,*5X 
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lUMBER OF TEETH OF SECDo GEAR=*»F12 0 6 ) 

38 F ORMAT ( / IOX 9 * T8E DIAMETERAL PITCH OF FIRST GEAR SET=* »F12„ 5 » // 
IXOXj* THE DIAMETERAL PITCH OF SECDo GEAR SET=*, F12 0 5 ) 

39 FORMAT ( / X0X>* HELIX ANGLE FOR THE FIRST GEAR SET=* , Fl2 0 5 , * dEG* 
1EES* s // 10XT* HELIX ANGLE FOR THE SECDo GEAR SET=*,F12.5,* DEGREE 
2S*) 

40 FORMAT ( / IOX 9 * THE TRANSMISSION R/.TIO IN THE FIRST REDUCT I0N=* » , 
lFl2»5»// 10X 9 * THE TRANSMISSION RATIO IN THE SECT » REDUCT ION=* , F55 
2 o 5 ) 

41 FORMAT ( / 10X,* T8E PERMISSIBLE ERROR IN THE MANUFACTURE=*=F12 a ' ) 

44 F ORMAT ( / 10XT * INITIAL VALUE OF R= #, F10 o 2) 

65 FORMAT! // 10X»* THESE ATE THE VALUES FOR COl* CR1 S CLI 5 CSl* 

1 CHI.* CT 1 9 CF1 FACTORS * ) 

66 FORMAT! / 10X 9 * THESE ARE THE VALUES FOR C02, CR2, CL2, CS2, 

1 CH2, CT2, CF2 FACTORS * ) 

67 FORMAT!// IOX 9 *THESE ARE THE MODIFIED OPTIMUM DESIGN VARIABLES** 

68 FORMAT! lOX* A7(1H-)} 

70 FORMAT! 8Fl0o5 ) 

71 FORMAT ( 7 F'6 0 2 ) 

72 FORMAT ! 1 3 ) 

73 FORMAT! / IOX 9 * THESE ARE THE VALUES FOR CF* SENSo FAVORS* ) 

74 FORMAT ( / 10X, 11F l0 o 5 } 

75 FORMAT! / lOX, * THESE ARE THE VALUES FOR THE DESIGN VARIABLES*) 

76 FORMAT! / 10X* * THESE ARE THE VALUES FOR THE MINIMUM VALUES 

1 FOR THE DESIGN Va7IABLES *) 

77 FORMAT! / 10X 9 * THESE ARE THE VALUES FOR THE SPAN OF THE DESIGN 
1 VARIABLES *) 

78 FORMAT! / lOX* * NO OF ITIRATI0NS= * 13 ) 

80 FORMAT!/ IOX 9 * THESE ARE THE VALUES FOR EPSl, EPS2, EPS3, EPS^*) 

91 FORMAT!/// 10X9*INPUT DATA *) 

92 FORMAT ( lOX » 10(lH-)) 

138 FORMAT!/ 10X,* THS MODULE FOR FIRST GEAR SET =*»F12 0 4,// 1 X,* TE 
IE MODULE FOR THE sSfcDo GEAR SET =*,F12 0 4 ) 

209 FORMAT! /// IOX 9 ^OPTIMUM GEAR DESIGN *) 

210 FORMAT! IOX 9 20(1H-) ) 

211 FORMAT! / 29X,* FIRST PINION *,2X»* FIRST GEAR *,4X »* SECDo PH*I 

ION *,4x»* SECDo GEAR */ ) 


212 


FORMAT! 

/ 

IOX 9 * NUMBER OF TEETH 

* 9 1 6 9 

9X916 , 

6 X 9 IX > 


1 

5X916 ) 






5x 9 FlO 

213 


FORMAT! 

/ 

IOX 9 *P C D IN M M« 


*»FlOo2, 

5X9 F 1 0 0 2 9 


1 O 

2 , 5X 

9 FIO 

0 2 ) 




5X ? 

214 


FORMAT! 

/ 

IOX 9 *FACE WIDTH IN 

M 

Mo * 9 F 1 0 0 2 9 

5x.FlCo2, 


1 

FlOo 2 , 

5X, 

FI 0 0 2 ) 



5x,F10o2, 

5X, 

215 


FORMAT! 

/ 

IOX 9 ^MODULES IN Mo 

Mo 

* 9 F 1 0 e 2 , 


1 F 1 0 0 2 » 

5X 9 

FlO 0 2 ) 




5X, 

216 


FORMAT! 

/ 

IOX 9 *HEL I X ANGLE I 

N DEC REE*.F10»2* 

5X*F10.2, 


1 

FlO 0 2 , 

5X, 

Fl 0 0 2 ) 



5X9 F 10 0 2 , 

5X 9 

217 


FORMAT! 

/ 

IOX 9 ^PROFILE CONTACT 

ratio*, F i 0 o 2 


1 

FlO 0 2 , 

5X9 

F 1 0 0 2 ) 



, 5 X 9 F 1 0 0 2 

9 5 X , 

218 


FORMAT! 

/ 

IOX 9 *FACE CONTACT 

RATIO *5 F 10c 2 


1 

FlO 0 2 9 

5X9 

Fl 0 « 2 ) 



5X9 F lQ 0 2 9 

5X 9 

219 


FORMAT! 

/ 

lOX, *P L VELOCITY 

IN 

M P S*9F10o2, 
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1 F 10 o 2 , 5 x, F 1 0 o 2 ) 

220 FORMAT ( / 10X* *H0RSE POWER IN W E /" R * ■> F12.2, 

1 F 12 o 2 , 5X» F 1 2 „ 2 ) 

221 FORMAT ( / 10X* *HORSE POWER IN 3 END I NG* , F 12 0 2 , 

1 F 1 2 o 2 , 5X, FI 2 e 2 ) 

222 FORMAT ( / 10X» *TRANSMl SSI ON- RATIO *, FlO. 2, 

1 F 10 o 2 j 5X* Fl 0 o 2 ) 

223 FORMAT ( / 10X * FALLOW MANF ERR IN M'M o **Fl0.4» 

1 FlO. 4 * 5 X» F 1 0 o A- ) 

224 FORMAT ( / 10X* *C0RREC FA '2 3 R *,FlO«2, 

1F10.2, 5x, F10 o 2 ) 

231 FORMAT ( X// lOX* ^GEOMETRY OF GEAR SET * ) 

232 FORMAT ( lOXT 20(1H-)) 

233 FORMAT!/ 10X* * BASE= *,Fl 0 .2, * [ M c *) 

234 FORMAT! / 10X, * HEIGHT =* , F10.2, * M M. *) 

235 FORMAT! / 10X»* AREA= *, FlO. 2,* SQ= METS *) 

241 FORMAT! /// 10X,* INSPECTION DATA * ) 

242 FORMAT! 10X, 16( 'lH- )) 

243 FORMAT!/// 10X» *FIRST PINION *, 5X,*FIRST GEAR *, 5X,* SECDo 

1 PINION *, 5X , * SECDo GEAR * ) 

300 FORMAT! //10X,*INPUT DESIGN DATA*) 

350 FORMAT ( IOXj 16! 1H-) ) 

301 FORMAT!/ 10X,*HORSE POWER=* » Fl 0 * 3 ) 

302 FORMAT!/ 10X,*TURBlNE SPEED 58 * »F8. 1,*RPM*) 

303 FORMAT!/ 10X ,*PROPELLER SPEED=* , F8 0 3 ,*RPM* ) 

304 FORMAT!/ 10X ^OVERALL REDUCTION RATIO=20,0*) 

500 FORMAT! //10X,* OBJECT FUNCTION * 5 F10.3 ,// , 10X»* PENALTY FUNCTION. 
1 *, FlO. 3 ,// ,10X»-* TOTAL FUNCTION *, F10.3 ) 

501 FORMAT! 'lOXr * DESIGN VARIABLES *, / * 3{10X, 6F16.6 / )) 

502 FORMAT! 10X, * CONSTRAINTS *,//, 10! 10X, 8F10.5 /)) 

503 FORMAT! 10X» * GRADIENT VECTOR*,//, 3( 10.X» 5F10.6 / )) 

READ 71* COl, CR1* CLl* CSl* CHl» CTl, CFl 

READ 71* C02, CR2 * CL2, CS2, CH2, CT2, CF2 
N= 11 
NC= 36 

PI= 4 o 0*AT AN ( 1,0 ) 

RAD= PI/ 180.0 

HP= 7 0 5*10 o 0**4 

PHI AN= 20« 0* RAD 

RPM1= 3000.0 

EG= 29.0*10.0**6 

C3= 1.67 

AK= 861.0 

GO= 1 o 1 8*lOC0**6 

SALB1= 2.25*10.0**5 

S ALB 2= 2.25*10.0**5 

UG= 0.20 

FATS= 1.5*10. 0**5 
SHEAR= 8.5*10. 0**4 


5x * 

F12.2= 

5X» , 

5X, 

Fl 2 o 2 = 

5X* 

5X * 

FlO .2 9 

5 X , 

5X* 

FlO o 4 , 

5 X » 3 


5X*F10.2, 5x, 
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zPi= ioo 0 o 
2P2= 75 o 0 
S 1 ADI- 3 QoO 
SIAD2= 3q o 0 
P ND 1= 3 0 8 5 
PND2= 2 o 16 
TRAN1- 5.0 
ALPHD= 65 o 0 
CNF 1= 11.5 
CNF2= 8 o 0 
ERR = 0 o 0008 5 
ALPHA= ALPHD*R AD 
SIA1= SIAD1* RAD 
SIA2= SIAD2* RAD 
Cl= C3*10.0**6*ERR 
TRAN 2= 20 o 0/ TRANl 

if if if if if if if if if -5c if if if “a- if if if if if if if if if if if if if if if if if ! if if 


* * **■ * ***#*•&#**#**#*•«•*■«■** * *** * 


THE DESIGN VARIABLES ARE NORMALISED AS UNDER 

if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if # if # if # if if if if if if if 

READ 70, CF, SENS 

READ 70, ( XMIN( I ) » 1= 1» N > 

READ 70, (Xl(I), 1= 1, N ) 

READ 70, ( SPAN ( I ) , 1= 1 , N ) 

READ 70, EPS 1 , EPS2, EPS3, EPS4 
READ 72, ITIR 


PRINT 

91 






PRINT 

92 






PRINT 

65 






PRINT 

74 , 

COl, 

CR 1 , 

CLl, CS1 , CHI, 

CT1 , 

CF 1 

PRINT 

66 






PRINT 

74, 

C02 , 

3 R2 , 

CL2, CS2 , CH2, 

CT2 , 

CF2 

PRINT 

73 






PRINT 

74, 

CF, 

SENS 




PRINT 

76 






PRINT 

74, 

( XM 

I N ( I ) 

, 1= 1* N ) 



PRINT 

75 






PRINT 

74, 

( xi 

( n , 

1 = • 1» N) ' A 



PRINT 

77 






PRINT 

74, 

( SP 

an( n 

, 1=1, N S 



PRINT 

80 






PRINT 

74, 

EPSl 

, EPS 

2,- EPS3 , EPS4- 



PRINT 

9\ 

00 

h- 

ITIR 






DO 505 1= 1, N 
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50- X 1 ( I ) = ( X 1 ( I ) - XMIN(I))/ SPAN (I) 

PRINT 79 

PRINT 74, ( XI ( I) » 1= 1 , N ) 

CALL PENAL! FO, FP » FT, XI, N» l o 0, GX ) 

DO 43 1= 1, NC 

43 I F ( GX(I)o LT o 0 o 2 ) FP = FP.+5 , 0-1 „ 0 / GX(I) 

PRINT 500, FO, FP, FT 
PRINT 502, ( GX ( I ) , 1= 1, NC) 

R= 0 o 06 17 
R= 5 o 0 

PRINT 44, R 
L3= „ 

I CALL UNC0NS ( XI* N, R ) 

DO 5 i= 1, N 

5 XN ( I ) = XI ( I) 

CALL PENAL( FO, FP , FT, XN, N, R, OX ) 

PRINT 500, FO, FP, FT 
PRINT 501, ( XN ( i ) , i= 1, N ) 

PRINT 502, ( GXCI) , 1= 1, NC ) 

CALL GRAD( XN, N, G, R ) 

PRINT 503, ( G ( I ) , 1= 1, N ) 

L 3 = L 3 + 1 ' 

I F { L3 0 EQo 1 ) GO TO 6 

I F ( ABSUFO- FLO ) / FO ) „ LT «, EPS4 ) GO TO 50 

9 SS = ABS( ( FT-FL) /FL) 

I F ( SSo GTo EPS1) GO TO 3 
DO 10 1= 1, N. 

A 1= ABS( ( XN ( I )- XL C I ) ) / XL ( I ) ) 

I F ( Alo GTo EPS2) GO TO 3 
10 CONTINUE ■ 

DO 11 1= 1, N 

II XO(I)= ( CF**0„5*XNt I)-XL(I).)/(CF**--0o5- l.o ) 
CALL PENAL! FO, FP , FTl , XO , N, R, GX ) 

I F ( FTlo GTo IOoQE+ 20) GO TO 3 
I F ( FUc GTo FT) GO TO 3 
SS= ABS( ( FTl- FT)/ FT) 

I F ( SSo GTo OoOOl ) GO TO 3 
DO 12 1= 1, N 

I F ( Alo GT, OoOOl ) GO TO 3 
Al= ABS(( X0(I)- X N ( T ) ) / XN(I)) 

12 CONTINUE 

GO TO 26 
6 RL= R 

DO 7 1= 1, N 

7 XL ( I ) = XN ( I ) 

FLO= FO 

FL= FT 
R= R/CF 
DO 8 1= 1, N 

8 XI ( I ) = XN ( I ) 

I F ( R 0 LEo 10 o 0E-3 ) GO TO 50 
GO TO l 
RL= R 


3 
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FU Ff 
Fl+O* fO 
ft* ft/CF 
00 15 I? 1, N 

xu n* xn( i ) 

i5 KB ill « KN( f 1*1*0/ CP*#6*5*(XN(H*XL(ti) 

FfttNt 16 j ( XR ( I ) ? 1= 1. N J 
CAfcU. P|NAL( f 0, FP» FT2 , XR? N,„R >(S GX ) 
IFVPT^# QT« 10o5E+2 n ) go TO 18 
I F ( ' PT^o GT. FL ) GO TO 18 
DO 17 J* 1, N 


19 

jam* xr ( 1 1 



I F ( R 6 LEb IObOEM ) GO TO 
GO TO lOQ 

50 

18 

DO 21 1= 1, 'N 

XL { I ) = XN( I ) 


21 

Xl ( I ) = XN( I ) 

GO TO 1 


50 

CONTINUE 

DO 51 1= 1, N 


51 

xot n= xi ( i > 


26 

PRINT 23, ( XO(I), 1= 1, N 
DO 24 1= 1, N 

5 

24 

XO(I)* XO( I )*SPAN( I )+ XMIN(I) 


PRINT 25, ( XO(I), 1= 1, N 
DO 30 1= 1, N 

) 

30 

XO( I ) = ( XO ( I ) - XM IN < I ) ) / 

SPAN ( n 


CALL PENALt FO» FP » FT? XO, N, R, GX ) 

PRINT 27, DPI? DGl 

PRINT 28, DP2, DG2 

PRINT 29, Bl, B2 

PRINT 31, PCR1, F3R1 

PRINT 32, PCR2, FCR2 

PRINT 33, PLVl, PLV2 

PRINT 3 A , HPWl, HPBl 

PRINT 35, HPW2 , HPB2 

PRINT 36, ZP1? ZP2 

PRINT 37, ZG1? ZG2 

PRINT 38, PNDl? PND2 

PRINT 39, SIAD1, SIAD2 

PRINT 40* TRAN1, T RAN 2 

PRINT 41, ERR 

CALL MODI F { XO, N, NC, R ) 

PRINT 67 
PRINT 68 

CALL DESIGN( XO, N, NC ) 

PRINT 27, DPI? DGl 
PRINT 28, DP2? DG2 
PRINT 29, Bl, B2 
PRINT 31, PCRl? FCR 1 
PRINT 32* PCR2, FCR2 
PRINT 33, PLVl, PLV2 


PRINT 
PRINT 
pjiUNT 
PRINT 
PpINt 
PRINT 
fP R I N T 
P'#INf 
PRINT 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
I 1 = 
ilp2= 
ttoi- 

I ZG2 = 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 

i-Rii'i r 

STOP 

END 


34-* NPWl. HPB1 
35, H P W 2 , HPB2 
3S» ZPl, £p 2 
* ZGl* ZG2 

iss » pmodI j pmodh 

39* SIAD1, SIAD2 

40, TRAN1 * Tran2 

41, ERR 
300 
350 

301, HP 

302, RPMl 
30-3, RPM3 
304 

209 

210 
211 
ZP1 
ZP2 
ZGl 
ZG2 

212, IZP1» IZG1» I ZP2 , izg2 
2139 DPI 9 DGI 9 DP2, DG2 

214, Bl, Bl* B2, B2 

215, PMODI 9 PMODl 9 PM0D2, PMC D2 
224, C0RRP1, CORRGI, C0RRP2 , C0RRG2 

216, Si ADI 9 SIAD1, SIAD2, SIAD2 

217, PCR 1 9 PCRl 9 PCR2, PCR2 
2189 FCR 1 9 FCR1 , FCR2 9 FCR2 

219, PLV1, PLVl , PLV2, PLV2 

220, HPW1 , HPW1 , HPW2 , HP W2 

221, HPB 1 , 8PB1, HPB2, HPB2 

222, TRAN 1 , TRANl, TRAN2, TRAN2 

223, ERR, ERR, ERR, ERR 

231 

232 

233, BASE 

234, HEIGHT 
akla 


SIBFTC CONST 

subroutine const 

DIMENSION XMINC15), XN(15), SPAN(.IS), Xl(15) 

COMMON/ DMS 1/ PI, RAD, NC 

COMMON/ DMS2/ COl, CRl* CH» CSl, CHl» CT1, CFl 
COMMON/ DMS3/ C02, CR2 , CL2, CS2, CH2, CT2 , CF2 
COMMON/ DMS4/ HP, PH I AN, EG, UG » RPMl, RPM2 , RPM3 
COMMON/ 0MS5 /DPI , DP2 , DG1 ♦ DG2 , B 1 , B2 , ZGl , ZG2 , PLVl , PLV2 , COSB 
COMMON/ DMS 6/ Cl, C3, AK, GO, CF, SENS, N 
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COMMON/ DMS7 / SALBl 
COMMQN/ DMS 8/ CDl, 
COMMON/ DMS9/ TORI, 
COMMON/DMSIO/ EPSl * 


COMMON/ DMS 11/ ZPl, 
COMMON/ DMS 12/ SIA1, 
COMMON/ DMS 13/ CNF1 
COMMON/ DMS 14/ XMIN 
COMMON/ DMS 16/ HPW1 , 
COMMON/DMS18/ CORRP1 
COMMON/DMS 19/ PMODl 
COMMON/ DMS 20 / FCR 
TRAN2= 20 o 0/ TRAN1 
DP1= ZPl/( PNDl*COs ( 
DP2= zP2/(PND2* COS ( 
B 1= CNF1*PI*C0S(SIA1 
B2= CNF 2* PI* COS ( S I 
DG1= DPI* TRAN1 
DG2= DP 2* TRAN2 
RPM2= PPM1/ TRAN 1 
RPM3= 1 50 o 0 


» SALB2, FATS, SHEAR 
CD2, CD^ » OZP 1 »OZGl » 0ZP2, 0ZG2 
T0R2, T0R3, BM1 , BM2 , B-3 
EPS2, EPS3, EPS4, ITIR 
ZP2, PNDl, PND2 , SIAD1, SIAD2 
SIA2, TRANl, TRAN2, ALPHD, ALPHA 

» CNF2 ’ ERR, DSH1, DSH2 , DSH3 , SCI, SC2 
XN, SPAN 

HPW2, I IP B 1 s, HPB2 , HEIGHT, BASE, AREA 
CORRG 1 , C0RRP2, C0RRG2 
TRAN,. FOB 

PCRl, PCR2 , WWl, 


, PM0D2, 
1» FCR2, 


WW2 , WD 1 , WD2 


SIA1) ) 
SIA2) ) 

)/ PNDl 
A2 ) / PND2 


TOR 1= 63000 o 0*HP/ RPMl 

TOR2= 63000. 0*HP/ RPM2 
TOR3 = 63000.0*HP/ RPM3 

BMl= 0 o 50* TORI 
BM2= oo 50* TOR2 
BM3= o o 50* TOR3 

******************** , ***********************-5f.#*i;-*-5{.*-jf.*******#j(.**#* 


CALCULATION of PITCH CONTACT RATIO AND FACE CONTACT RATIOS 

******************** *************** tt***#***#***#*#***.*^.^.^.^*^ 

PCRl = 0 o 555*TR AN 1**0 o 03 2* ZP 1**0 0 104* COS ( S I A 1 ) **1 0 5o4 
FCR 1 = CNF1* SIN(SIAl)* COS(SIAl) 

PCR2= 0.555* TRAN2**0 o 032* ZP 2**0 = 104* COS ( S I A2 ) **l . 5o4 
FCR2= CNF2*S IN ( S I A2 ) * C0S(SIA2) 

P LV1 = P I*DP1*RPM1/ 12.0 
PLV2= P I*DP2* RPM2 / 12 o 0 
WT1= 33000. 0*HP/ PL Vl 
WT2= 33000.0* HP/ PLV2 

IF(PLvl. GTo 5000o0)CVl- SORT ( 78 , 0/ ( 78 «. 0+SQRT ( PLVl ) ) ) 

I F ( PLVl. LE. 5000.0) CVl= 78.0/(78 .0+SQRT (PLVl ) ) 

I F ( PLV2 o GTo5000.0) CV 2= SORT ( 7 8 . 0/ ( 78 . 0+SQRT ( PL V2 ) ) ) 

I F ( PLV2.LE.5000.0 ) CV2= 78 . 0/ ( 78 . 0+SQRT ( PLV2 ) ) 

600 FORMATC // 10X* 3F15.5 ) 

BMAX 1= SORT ( 2 o 0*WT1 / ( ERR*GO ) ) 

BMAX2= SORT ( 2 o 0*WT 2 / ( ERR*GO ) ) 

I F ( BMAXl.GT.Bl ) CM1= BMAX1/ ( BMAXl-0 0 5*B1 ) 

I F ( BMAX2.LE.B2) CM2= 2.0*B2/ BMAX2 
I F ( BMAX2.GT.B2) CM2= BM AX2 / ( BMAX 2-0 . 5*b2 } 

GFl= 0.^5 o*TRAN1/(TRAN1+1,0) 

GF2= 0«^50*TRAN2/(TRAN2+1 o 0) 



n o 
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CP1= SQRT(EG/(2.0*PI*(1o0-UG*UG) ) ) 

CP2= CPI . 

SC1= CPl*SQRT (WT1*C01*CS1*CM1*CF1/ <CV1*B1*DP1*GF1) ) 

SC2= CP2* ,SQRT(WT2^C02*CS2'*CM2*CF2/(CV2*B2*DP2*GF2) ) 

SCDPl= ISALBl*DPl*CLl*CHl/(CPl*CTl*CRl) )**2 
SCDP2= (SALB2*DP2*CL2*CH2-/<CP2*CT2*CR2i )**2 
HPWl= RPMl*Bl*GFl*CVl*SCDPlO( 126 0 0( ,0*CSl*CMl*CFl*C01 ) 

HPW2= ( RPM2*B2/1260 a 0o 0 )* SCDP2/ ( CS2*CM2*CF2*C( .2 } * CV2 
WDT 1= ( B1*C1*C6S(SIA1 )*C0S(SIA1 )+WTl ) 

WDT2= ( B2*C 1*C0S { S I A2 ) *COS ( SI A2 ) + WT2 ) 

WD1= WT 1+ ( 0.« 05*PLV1*WDT 1*C0S ( Si A1 ) ) / ( SORT ( WDT 1 ) +0 o 05*PLVl> 

WD2= WT2+Oo 05#PLV2*WDT2*COS(.SIA2 ) / ( Oo 05*PLV2+SQRT ( WDT2 ) ) 

Ql= 2,0*TRANl/(TRANl+lo0) 

Q2= 2 o 0*TRAN2/ (TRAN2+1 o 0) 

WW1= PCR1*DP1*B1*A<*Q1/ (COS(SIAl)*COS(SIAl) ) 

WW2= PCR2*DP2*B2*AK*Q2/ ( COS ( SIA2 ) *COS ( SI A2 ) ) 

■YP1= 0.550-l o 41/(ZPl/(C0S(SlAl)**3)- l o 39)**0«756 
YP2= A0 55Ol 0 4-l/ (ZP2/ ( COS ( SI A2 ) **3 ) ~1 0 89 ) **0.756 

HPBl= ( RPM1**0 o 7 87*TRAN1**0 o 0-40/74000 o 0)*( FATS**1 „ 26/c3**0o067 ) ** 

1 ( DP1**3 o 224*cNF 1**1 o 9 2*YP1**1 o 26/ ( COS ( SlAl ) **2. 004* (ZP 1/COS ( SI Al 

2 ) **3 ) **2 O 450 ) ) 

HPB2= ( RPM2**0o787*TRAN2**0.0A0/74000,0)*( FATS**l a 26/C3**0o067)** 

1 (DP2**3.224*CNF2**1.192*YP2**1,26/ (COS ( SI A2 )**2 0 004*.( ZP2/COS (SI A2 

2 ) **3 ) ** 2 o 450 ) ) 

*****#*****M******************* ' * ***** ******* ****** ************** 

DSHl= < 16oO*SQRT(BMl**2+TORl**2)/(PI*SHEAR) )**0 o 333 
DSH2= ( 16 o 0*SQRT(BM2**2+T0R2**2 ) / (PI*SHEAR) )**0 o 333 
DSH3= ( 16 o 0*SQRT ( BM3**2+T0R3**2 ) / ( P I*SHEAR ) ) **0.333 
RETURN, 

END 


$ I BFT c fibo 

SUBROUTINE SEARCH ( STEP, XO, N» R ' 

COMMON/ DMS 17/ D* Y» G 

dimension xo(i5) • dc is > ,xm(15) »xki5) »X2<15> »gx(4oj »xtib) 


DIMENSION Y ( 15 ) , G( 15 ) 
INTEGER FN ( 12 ) 

data fn/i» i» 2 , 3 , 5 , 
DO lol I = 1 » N 
101 X ( I ) = X0( I) 

DIST= OolO 
1 Sub = DIST/21=0 

logical one, two 

ONE = o FALSE o 
TWO = o FALSE 0 
DO 20 I TER= 1 9 5 
NTEMP = 8 - ITER 
ND = FN(NTEMP) 


8, 13, 21, 34, 55, 89, 144/ 
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ALFA1 = SUB*FL0AT(ND) 

ALFA2 = SUB* FLO AT ( ND-1 ) 

DO 5 1 = 1, N 

XI ( I ) = XO ( I ) + ALFA1*D(I) 

5 X2 ( I ) = XO(I) + ALFA2*D ( I ) 

IF (ONE) GO TO 6 

cTall penal ( fp , fo, fi, xi, n. 

6 IF (TWO) GO TO 7 

call penalifo, fp, f2, x2, n, 

7 CONTINUE 

PRINT 900, ITER, Fl , F2 
IF ( Fl 0 LT 0 1 o OE+20 ) GO TO 9 
IF ( F2 „LTo lo OE+20 ) GO TO 8 
DIST = ALFA2 
GO TO 1 

8 DIST = ALFA1 

GO TO 1 

9 IF ( I TER o EQ o 5 ) GO TO 25 
900 FORMAT ( I 5 , 2E20.8) 

I F ( Fl 0 GT o F 2 ) GO TO 15 
DO 10 1 = 1, N 
10 XO ( I ) = X 2 { I ) 

F 2 = Fl 
TWO = 0 TRUEo 
ONE = « FALSE o 
GO TO 20 
15 Fl = F 2 

ONE = o TRUE o 
TWO = oFALSEo 
20 CONTINUE 

25 FMIN = AMIN1UF1» F2 ) 

I F < FM I N o EQ >. F2) GO TO 35 
DO 30 1 = 1, N 
30 XM( I ) = X1(I) 

GO TO 45 
35 DO 4o 1 = 1, N 
4-0 XM ( I ) = X 2 ( I ) 

45 step = (XM(i)-x(D )/dcd 
DO 50 1=1, N 
50 XO ( I ) = X ( I ) 

RETURN 

END 


R » GX ) 

R , GX ) 


SIBFTC PENAL 

subroutine penal< fo, fp, ft, x, n, r, gx j 
dimension x(i5), gx(4oj» xmincisj, xki5>, span(15>, xnci5) 
COMMON/ DMS 1/ PI, RAD, NC 

COMMON/ DMS2/ COl, CR1, CL1, CS1, CHI, CTl, CF1 
COMMON/ DMS 3/ C02, CR2 , CL 2 , CS2 , CH2 , CT2,CF2 


n n n n n 


31 


COMMON/ DMSA/ HP 9 PH I AN 9 EG 9 UG, RPM1 9 RPM2, RPM3 

COMMON/DMS5/DP1 *DP2,DG1,D+2,B1,B2,Z.G1 »ZG2,PLVl.PLV2,COSB 

'• i, “' rui '^ r7 ' SALBl» SALB2, FATS 9 SHEAR 

CD2, CD3 9 OZP 1 9 OZGl 9 OZP2, OZG2 
TOR2 ? TOR3 9 BMI 9 BM2, 0M3 
EPS2 , EPS3, EPSA, ITIR 
ZP2 , PNDl, PND2, SIADI 9 SIAD2 
* SIA2, TRAN 1 9 TRAN2, ALPHD 9 ALPHA 
9 CNF2 9 c. R R 9 D S H 1 9 DSH^9 DSNS 9 S C 1 9 SC2 
DMS 1 A/ XMIN, XN» SPAN 


C6MMON/DMS7/ 

C64M0N/ DMS 8/ CDl 9 
C6MMON/ DMS9/ TORI 9 
C6MMON/DMS10/ EPSl 9 
C6MMON/ DMS 11/ ZPl 9 
DMS 12/ S I A 1 
DMS 13/ CNF1 


C6MMON/ 

C6MMON/ 

C6MM0N/ 

C6MMON/ 

C6MMON/DMS18/ CORRPl, CORRG'l 
C6MMON/DMS 19/ PMOD1, PMOD2, 

COMMON/ DMS 20 / FCR1, FCR2 ? PCR1 9 PCR2 9 WW1, WW2 9 WD1» WD2 

* 1 *****************************M#********-* ****** ******** ********* 


DMS 16/ HPW1 9 HPW2 9 IIPB1, HPB2 


5 HE I GHT 9 
9 CORRG2 
TRAN 9 FOB 


■ 0 

CORRP2 


BASE, AREA 


CALCULATION OF CENTER DISTANCES ANl OBJECTIVE FUNCTION 


90 

91 


* 1 ************ ************** ****** ******** *********************** 


DO £ 0 1= 1, N 
XI ( I)= X( I) 

DO 91 1= 1, N 

X(I)= X(I)* SPAN (II*- XMl N ( + ) 

zpl= X(l) 

Z p2= X ( 2 ) 

PR D 1 = X(3) 

PKD2= X ( A ) 

SIAD1= X ( 5 ) 

SIA D2 = ^X(6) 

TR AN 1 = X ( 7 ) 
alphd= X ( 8 ) 

QJF1= X ( 9 ) 

CrF2 = X ( 10 ) 

EBR= X ( 11 ) 

CALL CONST 
TRAN 2= 20oO/ TRANl 
ZGl= OP 1* TRANl 
ZG2= OP 2* TRAN 2 


SENS= 0 o 50 

CD1= DPl # (TRA.Nl + lo0)/2 o 0 
CD2= DP2*(TRAN2+1 o 0)/2 o 0 

CD33= ( CDl**2 ) + ( CD 2** 2 ) 2 0 rt *CDl*Cf'2*CCS ( ALPHA ) 

CD3= SQRT(CD33) 

CO SB= ( ( CD1**2 ) + ( CD3**2 ) - ( CD2**2 ) ) / ( 2 „ 0*CDl*CD3 ) 
BETA= ARCOS(COSB) 

BET AR= BETA/ RAD 
BASE= CD3 

HEIGHT= CDl* SIN(BETAR ) 

AHEA= CD3*CD1* SIN (BETA )/(2 o 0* 1AA C 0 ) 

F0= ZPl*ZP2*{TRANl + l o 0)*( T-AN2+1 o _0 ) *Sl N ( ALPHA )/ ( 
1 1 )*COS(SlA2)*2 o 0) 


n 


PNDl*PND2*COS ( Slk 
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C 

c the following are the constraints for the problem 

c 

c 

G X ( 1 ) = 1.0- 1 3 0 7/ ZP1 

GX(2)= l.O- 13,7/ ZP2 

G X ( 3 ) = l.O-T.O/FCRl + 0.8*PCR1/ FCRl 

GX(A)= 1 oO-l o 0 / FCR2 + 0.8#'PCR2/FCR2 

GX(5)= 1.0-CNF1/ 12 o 0 

GX(6)= 1,0— A 0 0 / CNF 1 

G X ( 7 ) = 1,0- CNF2/ 12,0 

GX ( 8 ) = 1,0- A . 0 / CNF2 

GX ( 9 ) = 1,0- S I A D 1 / 35,0 

G X ( 1 0 ) = 1.0- 15.0/ SIAD1 

G X ( 11)= 1.0- 5 I AD2 / 35,0 

G X ( 12)= 1.0- 15,0/ SIAD2 

GX( 13)= 1,0- WDl/ WW1 

G X ( 1 A ) = 1.0- WD2 / WW2 

G X ( 1 5 ) = 1.0- HP/ 8 P W 1 

GX ( 16 ) = 1,0- HP/ HPW2 

GX( 17)= 1.0- H7 / HPB1 

G X ( 1 8 ) = 1,0- HP/ HPB2 

GX ( 19 ) = 1.0- (SC1*CT1*CR1)/ SALBl*CLl*CHl 
GX(2 a )= 1,0- (SC2#CT2*CR2)/ SALB2*CL2*CH2 
GX ( 2 1 ) = 1.0- ERR/ 0,0012 
GX ( 22 ) = 1.0- 0.0005/ ERR 
GX ( 23 ) = 1,0- 0 o 6*DG2/ CD3 

G X ( 2 A ) = 1.0- ABS ( 0 o 5*CD3~ CDl*COS (BETA) ) / ( 0,25* CD3 ) 
G X ( 2 5 ) = 1,0- B 1 / ( 6.0* PNDl) 

GX ( 26 ) = 1.0- (1.5*PND1)/ Bl 
GX ( 27 ) = 1.0- B 2 / ( 6.0*PND2) 

GX ( 28 ) = 1,0- ( 1 o 5*PND2 ) / B2 

GX ( 29 ) = 1.0- ( 2 . 0*DSH1 ) / DPI 

GX ( 30 ) = 1,0- ( 2 . 0"*DSH2 ) / DP2 

GX ( 3 1 ) = 1,0- ( 2 , 0*DSH3 ) / DG2 

GX ( 32 ) = 1.0- DP2 / ( 2 o 0 (CDl-15,0)) 

G X ( 3 3 ) = 1,0- ALPHD/90.0 
GX ( 3A) = 1.0- 30.0/ ALPHD 
GX ( 35 ) = 1,0- TRAN1/ 7.0 
GX ( 36 ) = 1.0- TRAN2/ 7,0 
X X = GX ( 1 ) 

DO 8 i= 1, NC 
I F ( XXoGE.GXd ) ) XX= GX ( I) 

8 CONTINUE 

I F ( XXoLT.0,0) GO TO 9 
SUM= n .O 
DO 10 1= 1, NC 
10 SUM= SUM+R/ GX(I) 

FP= SUM 
FO= FO**SENS 
FT= FO+FP 
GO TO 11 



9 FT= lOoOE+30 

11 DO 5 1= 1, N 

5 X ( I) = X 1 ( I) 

RETURN 
END 


;i bft c uncons 

subroutine uncons(X,n,r) 

dimension x c i 5 >■ , dusi, gx(40), y ( is > 

dimension gust, gkisi, h t is » 15), gxk40)» gx2(4o> 

COMMON/ DMS 10/ EPS 1 s EPS2 s EPS3 , EPS4, ITIR 
COMMON/ DMS 15/ GX , GXl, GX2 
COMMON/ DMS 17/ D* Y» G 

15 FORMAT ( // 10X, ^COUNTER NOT SUFFICIENT FOR MINIMISATION*) 

16 FORMATt // 1 OX s ^NORMAL CONVERGENCE AFTER* » I 3 ,* ITlRATlON*) 

25 FORMAT ( / lOX,* INITIAL DESIGN VECTOR * , 5 f 6 s 3 ) 

30 FORMATt// 1 OX , #M iNl M I SAT I ON OVER R *,Fl5 .6 , 3x »*STARTS*//J 

41 FORMATt// 10X, # LESS THAN THE LEAST VALUE OF X *) 

42 FORMATt / 10X, * GRADIENT =*,6El5„6 ) 

43 FORMATt / 10X»* DIRECTION =* 6E15„6 ) 

44 FORMATt 10X,* STEP-SIZE ■* » E15o6 ) 

49 FORMATt / 10X»* FROM UNCONS #, 5X, * FO =*,E15 0 6,* FP =*,E15 0 6, 

1 * FT1 = *, El5 o 6 ) 

50 FORMATt 10X, * ANGLE = #, F6 0 3 ) 

51 FORMATt/ 10X, * MODIFYING MATRIX * ) • 

132 FORMAT (15X»*DERIVATIVE OF FUNCTION ALONG DIRECTION OF *,/»15x, 

1*LINEAR MINIMISATION IS NOT NEGATIVE AT*, 14,* IT I RAT ION FOR R = *-»0 

2 F 1 0 o 6 ) 

PRINT 30 ,R 

PRINT 25 , (X( I ) , I=1»N) 

L7 = 0 
LI = n 

CALL PENAL(FO,FP,FT,X»N,R,GX) 

I F ( FT o GT o lOoOE+20 ) GO TO 17 

call graD( x,n,g,R) 

PRINT 42, (G( I ) *I=1»N) 

1 DO 5 13 = 1 ,N 
DO 5 12 = 1 , N 
H< 12, 13) = OoO 

5 H( 13, 13) = loO 

8 DO 6 II = 1,N 

D( U) = 0 o 0 
DO 6 12 = 1,N 

6 D ( 1 1 ) = D<I1)-H(I1»I2)*G(I2) 

S9 = oo 0 

S8 = 0 o 0 
DO 10 I = 1»N 
S 8= S8+D(I)*D(I> 

10 S9 = S9 + D ( I ) *G t I ) 
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IF(S9 o LTo0o 0) GO TO 2 
DO 23 1 2 = 1 , N 
I F ( H ( I 2 » I 2 ) o NE o 1 o 0 ) GO TO 1 
23 CONTINUE 

PRINT 132, LI, R 
RETURN 

2 XP = l o 0 

DO 199 I = 1 , N 
I F C X C I ) .LToXP) XP = X ( I ) 

199 DC I ) = D ( I ) /SORT (S8 ) 

PRINT 43, ( n(i), i= 1, Nj 
XP = XP*0o00l 

call search(Sig,x,n,r) 

IF( SlGo LEo 10 o 0E-6 ) R0 TO 13 
SlG = SIG/SQRT(S8) 

PRINT 44,SIG 
IF CSIGoLT„XP) L7 = L7+1 
I F ( L7o GTc 10 ) GO TO 40 
DO 7 11= 1 , N 
D( II) = DC II )*SQRT(S8) 

X ( 1 1 ) = X ( 1 1 ) +S I G*D ( I 1 ) 

7 Gl(Il) = GCIl) 

PRINT 25, (XCI) , 1 = 1, N) 

LI = Ll+1 

IF (Ll„GTolOO) GO TO 14 

call penal t fo,fp,fti ,x»n,r,gx) 

PRINT 49,F0,FP,FT1 
Ss = ABS( (FT1-FT) /FT) 

FT = FT 1 

call grad , cx»n,g»r) 

PRINT 42, (G( I ) , 1=1, N) 

ET = AMAX1 ( loOsSIG) 

L2 = n 

DO 3 I = 1,N 

IF(ABS(ET*D( I ) ) oLToOoOOl) GO TO 3 
L2 = L2+1 

3 CONTINUE 
IF(L2 0 FQoO) GO TO 13 


SI 


A 0 0 

S2 


A « 0 

S3 

= 

A 0 0 

DO 

11 

I = 1 , N 

si 

= 

Sl+GC I 

S2 

= 

S2 + D ( I )* *2 

S3 

= 

S3+G( I )*D( I ) 

S2 

= 

SOFT ($2 ) 

SI 

= 

SOFT (Si ) 

SG 

= 

ARS ( S3/S2/S1 ) 


PRINT 50, SG 

IF (ABS(S3/Sl/S2) o GTo0oD GO TO 8 
cxcx = loOE-2 

IF (S1oLEoCXCXoOR 0 S2oLEoCXCX) GO TO 13 



o n o o n o 
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151 

DO 9 ] 

tl = 1, 

N 


9 

Y{ 11 ) 

= G ( 1 1 

i — f 

0 

1 

u 


Sll = 

Do 0 




S12 = 

Go 0 




DO 15< 

-1 = 1 

, N 



Sll = 

S 1 1 + D ( 

I )*Y( 

i ) 

152 

S12 = 

S12+Y ( 

I )*G( 

i ) 


sum 2 = 

= OoO 




DO 156 I = 1 

, N 



SUM1 = 

= OoO 




DO 157 J = 1 

9 N 


157 

SUM1 

= SUM1 

4* X ( J ) 

*H 


156 SUM2 = SUM2+X ( I )*SUM1 
A = SUM 2 

IF(AoLt o 0 o 0) GO TO 1 
PRINT 51 

call hmatrt <h,sig,N) 

60 TO 8 

1 A PRINT 15 

GO TO 17 

40 PRINT 41 

GO TO 17 

13 PRINT 16, LI 

17 RETURN 

END 


SIBFTC MODIF 

subroutine modifc xo, n» nc, r ) 

DIMENSION XO (15), XN<15), SPANU5M X(15) 

DIMENSION P NDA (30), PNDB(30), XMIM15) 

DIMENSION GX ( 40 ) , 6X1(40), GX2(40) 

C0MM0N/DMS5/DP1 ,DP2,DGl ,DG2,B1,B2,I61,2 g2,PLV1,PLV2,cOSB 
COMMON/ DMS 8/ CDl, CD2, CD3 ,0ZP1 ,OZGl » OZP2,.OZG2 
C0MM0N/DMS11/ ZP1, Zp2 , PNDl, PND2, S'lADl, SIAD2 
COMMON/ DMS 12 / Sill, SIA2, TRAN 1 , TRAN2 , ALPHD , ALPHA 
COMMON/ DMS 13/ CMFl, CNF2, ERR, DSHl, DSH2, DSH3, SCI, SC2 
COMMON/ DMS 14/ XMlN, XN, SPAN 
COMMON/ DMS' 15/ GX , GX1..GX2 
COMMON/DMS 19/ PM0D1, PM0D2 , TRAN, FOB 

**************■«■***«•* ***-»**'-)t#*’#*Q***'********** ******************** 

the following are the optimum design variables without rounding 

OFF the NO of TEETH AND WITHOUT USING THE PREFERRED MODULES 

******** *****•}(•** ******%****■){•* **^ ******* ******^**** ************ 

DO 90 i= 1, N 
90 XO(l)= XO ( I ) * SPAN ( I ) + XMlN(I) 

ZP1= X0(1) 

ZP2= xO(2) 



U U U U U U U KJ KJ KJ 
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10 


PND 1= XO ( 3 ) 

PND2= XO ( 4 ) 

S I A D 1 = XO ( 5 ) 

S I A D 2 = XO ( 6 ) 

TRAN 1= XO ( 7 ) 

A L P H D = XO ( 8 ) 

CNF 1 = XO ( 9 ) 

CNF2= X0(10) 

ERR= XO(ll) 

PMODl= 2 5 0 4/ PND1 
PMOD2= 25 0 4/ PND?. 

TRAN2= 20 o 0 / TRAN1 
ZG 1= ZPl* TRAN1 
ZG2= 7 P 2* TRaN 2 
0ZP1= Z.P1 
0ZP2= ZP2 
0ZG1= ZG1 
0 Z G 2 = ZG2 

READ 10, ( PNDA ( I) , 1 = 1 , 30 ) 
READ 10, ( PNDB ( I ) , 1 = 1 5 30 ) 
FORMAT* 8 F lOo 5 ) 

ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft 


# -X” * ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft ft 


ROUNDING OF THE NUMBER OF TEETH OF PINION AND GEAR 


* **#■)£■*•!<■*■*(- #tt * ** **■»•** ********* ****** -k *****##*■*#*•£*#** -Jf. *##■* # * x 

IZP1= ZP1 
AZP1= IZP1 

I F ( ( ZPl-AZPl ) o GEo 0 o 5 ) AZPl= AZPl+ 1„0 

I ZP2= ZP2 
AZP2= I ZP 2 

I F ( ( ZP2-AZP2) 0 GE 0 0 o 5 ) AZP2= AZP2+ l„o 
IZGl= ZG 1 
AZGl= I ZG 1 

I F ( ( ZGl-AZGl ) o GEo 0 o 5 ) A Z G 1 = A Z G 1 + loO 
I Z G 2 = ZG2 ' 


AZG2= IZG2 

I F ( ( ZG2-aZG2) 0 GEo 5 ) AZG2= AZG2+ 1.0 


RE- STORING THE ROUNDED OFF VARIABLES IN THE ORIGINAL NAME 

#*-)£■•«•*• jh*- •*##*# ***##*• i tt***##********#*-#**##*##*##**-*#*#*****##*****# 

ZP1= AZP1 

ZP2= AZP2 

ZG1= AZG1 

ZG2= AZG2 

TRAN 1= ZG 1 / ZP1 

TRAN2= ZG2 / ZP2 

TRAN= TRAN1* TRAN2 
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15 

16 


17 

18 


20 


24 

30 


26 

35 


27 

4-0 


XO ( 1 ) = ZP1 
XO ( 2 ) = ZP2 
XO ( 7 ) = TRAN1 
DO 15 1= 1, 30 
J= 1+ 1 
11= I 
Jl= 11+ 1 

IF(PMODl<> LEo PNDA( J) 8 ANDo PMODl 0 GEo P N DA ( I } ) GO TO 16 

CONTINUE 

PMOD 1= PNDA(J) 

DO 17 1= 1, 30 
J = I + 1 


I 2= I 
J 2 = 12+ 1 


I F ( PMOD2c LEo P N D B ( J ) o AND 0 PM0D2, GEo PNDB(I)) GO TO 18 

CONTINUE 

PMOD2= PNDB(J) 

PND1= 25 0 4/ PM0D1 
PND2= 2 5 „ 4/ PMOD2 


XO C 3 ) = PND1 
XO ( 4 ) = PND2 
DO 2 n L= 1, N 

XO ( L ) = ( X0(L)-XMIN(L) ) / SPAN(L) 

X ( L ) = XO(L) 

CALL PENAL ( FO , FP» FT, X, N» R, GX ) 

DO 24 K= 1, NC 

I F ( GX(K) 0 LToOoOo ORo GX(K)o GTo l o 0 ) GO TO 3p 

GO TO 25 

PMOD 1= PNDA ( 1 1 ) 

PMOD2= PNDB ( 1 2 ) 

X ( 3 ) = 2 5 o 4/ PMODl 
X ( 4 ) = 2 5 „ 4/ PMOD2 
X ( 3 ) = ( X ( 3 ) — XMlN ( 3 ) ) / SP AN ( 3 ) 

X ( 4 ) = ( X ( 4 ) - XM.lN ( 4 ) ) / SPAN ( 4 ) 

CALL PENAL ( FO , FP , FT, X, N, R, GX ) 

DO 26 K= 1, NC 

I F ( GX(K)o LTo 0 o 0 o ORo GX ( K ) = GTo loO ) GO TO 35 

GO TO 25 

PMOD 1 = PNDA ( I 1 ) 

PMOD2= PNDB ( J2 ) 


X ( 3 ) = 25 0 4/ PMODl 
X ( 4 ) = 25 0 4/ PMOD2 
X ( 3 ) = ( X ( 3 ) — X M I N ( 3 ) ) / SPAN ( 3 ) 

X ( 4 ) = ( X ( 4 ) *- X M I N ( 4 ) ) / SPAN (4) 

CALL PENAL ( FO» FP, FT, X, N , n R , n GX n ) 

DO 27 K= 1, NC 

I F ( GX(K)o LTo OoOo ORo GX(K)o GTo loO ) GO TO 40 

GO TO 25 

PMOD 1= PNDA(jl) 


28 

25 

55 


76 


75 

60 


PMOD2= PNDB ( 1 2 ) 

X(3)= 25 0 4/ PM0D1 

Xf4) = 25 0 4/ PM0D2 

XO)= ( X ( 3 ) - XM I N ( 3 ) ) / 

X(4)= ( X ( 4 ) — X M I N ( 4 ) ) / 

CALL PENAL ( FO , FP , FT, 

DO 28 K= 1, NC 

I F ( GX(K ) 0 LT 0 OcOo OR 0 

GO TO 25 

CONTINUE 

GO TO 75 

11= J1 

Jl= J 1 + 1 

12= J2 

J2= J 2 + 1 

I F ( Jl 0 GE 0 28 0 OR „ J2. 

GO TO 30 

CONTINUE 

I ZG2= ZG2 

AZG2= IZG2 

I F ( ( ZG2-AZF2 ) 0 GE o 0 o 5) 
ZG2= AZG2 
DO 60 1 = 1 , n 
XO( I ) = X { I ) 

DP1= ZPl* PMOD1 
DP2= ZP2# PMOD2 
DGl= ZGl* PMOD1 
DG2= ZG2*PMOD2 
RETURN 
END 


SPAN (3) 

SPAN ( 4 ) 

X s N » R , GX ) 

GX(K)o GT o loO ) GO TO 55 


GE« 28 ) go TO 76 


AZG2= AZG2*p1 o 0 


SI BFT C DESIGN 

SUBROUTINE DES 
DIMENSION XO( 
COMMON/DMS5/DP 
COMMON/ DMS 8 / 
COMMON/ DMS 11/ 
COMMON/. DMS 12/ 
COMMON/ DMS 13 
COMMON/ DMS 14 
COMMON/ DMS 16 / 
COMMON/DMS18/ 
COMMON/DMS 19/ 


, XT 15), XN ( 15 ) 

9 ZG2 , PLVl »PLV2 ,C0SB 


DO 9 0 l= 1, n 
90 XO( I )= XO ( I )# 
ZP1= XO(l) 
ZP2= X0(2) 
PND1= XO ( 3 ) 

P N D 2 = XO ( 4 ) 


I GN ( XO, N, NC ) 

15 ) » XMIN{ 15) , SPAN( 15) 

1 >DP2,DGl ,DG2 ,B1,B2 ,zg 1 
CDl , CD2, CD3 s OZP 1 , OZGl » OZP2, OZG2 
ZPl, Z 72 , PNDl , PND2, SIADl, 5IAD2 
SlAl, SIA.2, TRANl, TRAN2, ALPHD, ALPHA 
/ CNF 1 , CNF2, ERR, DSH1, DSH2 , DSH3. SCI, 

/ X M I N , XN, SPAN ' 

HPW1 , HPW2 , I IP 8 1 , HPB2 9 HEIGHT, BASE, AREA 
CORRP1, C0RRG1, CORRP2 , CORRG2 
PMODl, PMOD2 , TRAN, FOB 


SC2 


SPAN ( I ) + XMIN ( I 


n n n n n 
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SIADl= XO ( 5 j 
S I AD2= XO ( 6 ) 

TRANl= XO ( 7) 

ALPHD= XO(8) 

CNF1- XO ( 9 ) 

CNF2= XO(10 ) 

ERR= XO (11) 

ALPHl= ALPHD*3 o 1416/ 18o a O 
ALFA= 20 o 0* 3 e 1416/ 180 o 0 

******* * * -5!- * ** -X- ** * ** * ** -if- * ** * * * ** * * * * * *** * * -ft * * * * * ** * ****.*■# * * * * *** -ft * 


ROUNDING OFF THE CENTRE DISTANCES 


C 

C 

C 

C 

C 


15 


******** * * * * * * * * * * * * * * * * 

CEN1= CD1 
CEN2= CD2 


* * ************ * * * * * * * *************** * * * * 


I CEN 1= CEN1 
ACEN 1= I CEN 1 

I F ( (CENl-ACENl ) « GE 0 0„5 ) ACENl= ACEN1+ 1.0 
ICEN2= CEN2 
ACEN2= ICEN2 

IF ( ( CEN2-ACEN2 ) 0 GEo 0.5 ) ACEN2= ACEN2+ 1.0 

CEN1= ACE.N1 

CEN2= ACEN2 

CD1= CEN1 

CD2= CEN2 

CD3 3 = (CD 1**2 ) + ( CD2**2)-2 o 0*CDl*CD2*C0S(ALPHA) 

CD3= SORT ( CD3 3 ) 

COSB= ( (CD1**2 )+(CD3**2)-(CD2**2) ) / ( 2 o 0*CDl*CD3 ) 

BET A= ARCOS(COSB) 

BASE= CD3 

HEIGHT= CDl.tF SIN(BETA) 

AREA= CD3*CDl*SIN(BETA)/(2e0*lA4 o 0 ) 

FOB= AREA 

********************************* ******************* ****** ******* 


ESTIMATION OF TOTAL CORRECTION, XT= 2x 

***************************************************************** 
Q 1= (0ZP1+ OZGl ) *COS (ALFA) / (ZPl+ZGl ) 

A31= aPCOS(QI) 

FORMAT ( / 10X> 2 El5 0 4 ) 

PRINT 15, Ql, A31 

XTl = ( TAN ( A3 1)-A3 1-TAN ( ALFAJ+ALFA )*(ZPl+ZGl )/ (2.0*7 AN (ALFA) ) 

Q2= (OZP2+OZG2 )*COS( ALFA) / (ZP2+ZG2 ) 

A32= APCOS ( Q2 ) 

PRINT 15, Q2, A32 

XT2= ( TAN ( A3 2 ) -A32-T AN ( ALFA ) +AL FA ) * ( ZP2+ZG2 ) / ( 2 . 0*TAN ( ALFA ) ) 



nnn 


‘90 

CgU p CpRREC< IPl ? ZGl? JRANi* XTl, XP, XG ) 

CCRRG1= XG 

CALL dORRECt t?$ 9 Zg2 » TRA.N2, XT2» XP , XG ) 

C0RRP2= XP 
CQPPG2= XG 
Vl- PLV 1 / 3 , 28 
V2= PLV2/3 0 28 

xi- coRrpi 

X2= CORRGl 
%%= C0RRP2 
X4= C0RRG2 

cOrrpi- xi 

C0RRG1= X 2 
C0RRP2- X 3 
C0RRG2= X4 
I QLTY= ERR 
I QLTY= 7 

*********** ** *** * ****** * **** * *** * **** * ** -x- ** * * * * * -X-#* * ** * *** ***** * 

CALCULATION OF MEASUREMENT DATA 

****** -if- ** * * * -if *** ** **** -St- * * * -X- * * * * -if- ** * * ***** ** * ** * ************ * **** * 

DPl = 25 o A* DPI 

DP2« 25.4*DP2 

DGl* 25 0 A* DGl 

DG2= 25,4* DG2 

B 1= 81* 25,4 

B2= 82* 25,4 

PLV1= PLV1/ (3.28*60,0) 

PLV2= PLV2 / ( 3 o 28*6C ® 0 ) 

ERR= ERR* 25 0 4 
BASE= BASE* 25,4 
HEIGHT* HEIGHT* 25,4 
AREA* AREA/ ( 3 o 28*3 ,28 ) 

RETURN 

END 


$IBFTc CORREC 

SUBROUTINE CORREC ( ZP , ZG* TRNS » XT , XP , XG ) 

DIMENSION XM ( 12 ) ? SLP ( 12 ) > XM1(13), SLPK13) 

COMMON/ DMS 16/ HPW1. HPW2, IIPB1, HPB2, HEIGHT, BASE, AREA 
DATA XM/~0 « 35 , -0.25, -0.14, -0.05, .0.08, 0.17, 0 = 26 , 0.38= 0.50, 
1 0.62, 0 o 72 , 0 o 84 / 

DATA SLP/-0.8935 , -0.7535, -0.6o85, -Q.A835, -o„ 3775, -0.2735, 

1 -0 o 1762 , -0.0906 ,-0o00A36, +0.0567, +0.01351, +0.1942 / 

DATA XMl/ -0 o 38, “0.18, -0.12, - 0 . 06 , 0.00, 0.07, 0.15, 0. 23 , 

1 0.34, 0. A4, 0.57, 0.72, 0.?. rt / ^ . 

DATA SLPI/ -0. 3 0 6 0 5 ~0.0 559 , ~0.oA16, -o«ol 7,+ » 0»0» 0.00 332 , 
10.0663, 0 o 1 0 5 0 * 0.162Q, 0.2325 , o* 33 0 3 , O.A75 0 , 0.685 0 / 



X = XT/ 2 o 0 


Si 


C 

C 

C 

C 

C DISTRIBUTION OF CORRECTION 

ZZ = ( ZP + ZG)/ 2 o 0 
,IF( TRNSo GT 0 1,0 ) GO TO 731 

c for speeding down pair of gears 

XX1= SLP (■ 1 ) * ( ZZ-40 0 0 ) / 50,0+ XM(1)- X 
DO 735 N= 2, 12 

XX2= SLP ( N ) * ( ZZ-40, 0)/ 50,0+ XM(N)- X 
I F ( XX1+XX2 ) 736, 736, 734 

734 XX1= XX2 

735 CONTINUE 

736 SL- S[_P ( N — 1 ) —X X 1 * ( SLP ( N ) —SLP ( N— 1 ) )/(XX2— XXll 

GO TO 751 ' 

C FOR SPEEDING UP PAIR OF GPARS 

731 IFfX.GEo (-0,3) V GO TO 732 

I F ( ZZ o GT, 5o,0 ) GO TO 739 

732 XX1= SLP 1 ( 1 ) * ( ZZ-40,0)/ 50,0+ XMl(l)- X 
DO 795 N= 2, 13 

XX2= SLP1 ( N )#( ZZ-40, 0 ) / 50,0+ XMKN)- X 
IF(XX1* XX2 ) 737, 737, 738 

738 XX 1= XX2 
795 CONTINUE 

737 SL= SLP 1 ( N— 1 ) “XXI* ( SLP 1 ( N ) - SLPl(N-l))/ ( XX2-XX1 ) 
751 XP= SUHZP-ZZ) / 50,0+ X 

XG= XT- XP 
GO TO 761 

739 XP= X 
XG= X 

761 RETURN 
END 


SIBFTC GRAD 

subroutine grad( x, n, g, p.) 

dimension x(i5j, cxi(40), gx2(4oi, cn.5) 

dimension gx(4o) 

COMMON/ DMS 1 / PI, RAD, NC 
COMMON/ DMS 15/ GX , GXl, GX2 

call penal r fo, fp, ft, x, n, r, g>U j 
Do 2 1 = 1 , n 

A= X < I ) 

ST= AMiN1( 0,001, A/10,,0,0 ) 

X( I )= X( I )+ ST 

call penal ( foi, fp, ft, x, n, r, gx2 > 

G( I )= ( FOl-FO) / ST 
DO 3 J= 1, NC 


G(I>= GUI- (( GX2 ( J ) — GX1 ( J) ) /ST) /GX1 ( J)**2 *R 
3 CONTINUE 

X ( I ) = X C I ) — ST 
2 CONTINUE 

RETURN 
END 


$ I bft c hmat 

subroutine hMatri < h, sig. n ) 

COMMON/ DMS 17/ 0, Y> G 

dimension hsiisi, H9(i5, 15 j » hioxjs, is), hshs) 
dimension H6 ( 1 5 » 15 ) > H(15,15), Y(lf), DC 15 } , HAPS) 
DIMENSION G(15) 

c 

Hl= a sO 
DO 1 1= 1, N 
HI = H1+ D( I)* Y( I ) 

H5 ( I ) = 0 o 0 
HA ( I ) = 0 o 0 
DO 1 J= 1, N 

HA( I ) = HA ( I D+H ( I »J ) *Y.( J ) 

H5( I )= H 5 ( I) + H( J,I)* Y( J) 

1 H9 ( I » J ) = D ( I ) * DC J ) 

HI 1= 0 o 0 

DO 2 1= 1, N 
H 1 1= H11+ Y( I )* H A ( I ) 

DO 2 J= 1, N 

2 HIO ( I , j) = HA< I )# H5 ( J ) 

DO 3 1= 1, N 

DO 3 J= 1, N 

H6(I» J)= H ( I » J ) + SIG/ H1*H9(I,JJ- HlO(I,J)/ Hll 

3 CONTINUE 

DO 5 1= 1, N 

DO 5 J= 1, N 

5 H( I, J ) = H6 ( I j J) 

RETURN 

END 


SENTRY 



SAMPLE' OUTPUT FOR EXAMPLE -2 


INPUT. DATA 

values for coi, cri 5 cli, csi 5 chi, cti ? cfi factors 

1.5000 1.2500 1.1000 1.0000 1.0000 1.0000 1.0000 

VALUES FOR C02 , CR2, CL2, CS2, CE2 , ,,CT2 , CF2 FACTORS 
1.5000 1.2500 1.1000 1.0000 1,0000 1.0000 1.0000 

VALUES FOR CF , SENS. FACTORS 

10.0000 1.0050 

VALUES FOR THE MINIMUM VALUES FOR THE DESIGN VARIABLES 
13.7000 13.7000 0.5000 0.5000 15,0000 15.0000 2.0000 

30.0000 5.0000 5.0000 0,0005 

VALUES FOR THE DESIGN VARIABLES 


100.0000 

130.0000 

2.0000 

2.0000 

28,0000 

28.0000 

5,0000 

70.0000 

8.0000 

8.0000 

0.0008 


VALUES FOR 

THE SPAN 

OF THE DESIGN VARIABLES 


250.0000 

250.0000 

15.0000 

15.0000 

20.0000 

20.0000 

5.0000 

60.0000 

10.0000 

10.0000 

0.0007 



VALUES FOR EPSl, EPS 2, EPS3, EPS4 
0;0005 • 0.0005 0.0005 0.0005 

NO OF ITERATIONS = 200 
OBJECT FUNCTION 242.714 
PENALTY FUNCTION 81.547 
TOTAL FUNCTION 344.492 



CONSTRAINTS 


0.86300 0.80462 0.88680 0.8S037 0.33333 0.50000 0.33333 0.50000 

0.20000 0.46429 0.20000 0.46429 0.90473 0.77785 1.00000 0.81678 

0.96397 0.52811 1.00000 0.41812 0,29167 0.41176 0.07274 0.72398 

0.09310 0.72434 0.09310 0.72434 0,84827 0.80053 0.92087 0.76279 

0.22222 0.57143 0.28571 0.42857 

INITIAL VAILS OF H = 200.00 
MINIMISATION OVER I = 2.00.00000 STARTS 
NORMAL CONVERGENCE AFTER 25 ITERATION 


OBJECT FUNCTION 590.172 
PENALTY FUNCTION 14305. 574- 
TOTAL FUNCTION 14895.746 
MINIMISATION OVER R = 20.00000 STARTS 
NORMAL CONVERGENCE AFTER 23 ITERATION 


OBJECT FUNCTION 251.449 
PENALTY FUNCTION 1633.859 
TOT AL FUNCTION 1786.308 
MINIMISATION OVER R = 2.00000 STAGS 
LESS THAN THE LEAST VALUE OF X 


OBJECT FUNCTION 150.152 ■ 

PENALTY FUNCTION' 187.348 

TOTAL FUNCTION 337.500 

MINIMISATION OVER R = 0.20000 STARTS 



NORMAL CONVERGENCE AFTER 20 ITERATION 

OBJECT FUNCTION 109.515 
PENALTY FUNCTION 31.404 
TOTAL FUNCTION 140.919 
MINIMISATION OVER R = 0.02000 0 STARTS 
NORMAL CONVERGENCE AFTER 15 ITERATION 

OBJECT FUNCTION 98.851 
PENALTY FUNCTION 6.405 
TOTAL FU MOTION 105.265 
MINIMISATION OVER R = 0.002000 ST .'UTS 
LESS THAN' THE LEAST VALUE OF X 

OBJECT FUNCTION 95.882 

PENALTY FUNCTION 1,572 

TOTAL FUNCTION 97.453 

MINIMISATION OVER R = 0.000200 STARTS 

LESS THAN THE LEAST VALUE OF X 

OBJECT FUNCTION 95.640 
PENALTY FUNCTION 0.218 
TOTAL FIJI CT ION 95.859 

NORMALISED OPIIuUk DESIGN VARIABLES 
0.087061 0.254675 0.123583 0.122399 0.399731 
0.399727 0.801157 0.373471 0.375708 0.687293 
0.012267 



96 


ACTUAL DESIGN VARIABLES 

35.4652 77.3688 2.3537 2.3360 22.9946 
22.9945 6.0058 52.4083 11.7571 11.8720 
0 .0005 


INPUT DESIGN DATA 
HORSE POWER = 60000.000 
TURBINE SPEED = 3000 .ORRi 
PROPELLER SPEED = 150.000RPM 
OVERALL REDUCTION RATIO =20.0 


I 1 1. T. KANPUR, 

j CENTRAL LIBRARY# 



OPTIMUM GEAR DESIGN 


NUMBER OF TEETH 
P C D IN M M. 

FACE WIDTH IN MM. 
MODULES IN M.M. 

CORREC FACTOR 
HELIX ANGLE IN DEGREE 
PROFILE CONTACT RATIO 
FACE CONTACT RATIO 
P L VELOCITY IN M P S 
HORSE POWER IN WEAR 
HORSE POWER IN BENDING 
TRANSMISSION RATIO 
ALLOW MANF ERR IN M M. 


FIRST 

FIRST 

SECOND 

SECOND 

PINION 

GEAR 

PINION 

GEAR 

34 

212 

76 

254 

442.00 

2756 .00 

912.00 

3048.00 

415 .84 

415.84 

387.63 

387.63 

13.00 

13.00 

12.00 

12.00 

-0.072 

5 -.158 

-.42 

-.42 

22.99 

22.99 

22.99 

22.99 

1 t 

0.69 

0.69 

0.73 

0.73 

.5.09 

5.09 

5.14 

5.14 

82 . 55 

82 .55 

27 .55 

27.55 

128208.55 

128208.55 

320922.92 320922.92 

1521862.27 

1521862.27 

93157.83 

93157.83 

6.09 

6.09 

3.29 

3.29 


0.0129 0.0129 


0.0129 0.0129 



GEOMETRY OF GEAR SET 


BASE = 
HEIGHT= 


1868.79 MM 
1797.23 MM 


AREA 


1.68 SQ.METS 



